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Abstract 

The Nakanishi perturbative integral representation of the four-dimensional T-matrix is investi- 
gated in order to get a workable treatment for scattering states, solutions of the inhomogeneous 
Bethe-Salpeter Equation, in Minkowski space. The projection onto the null-plane of the four- 
dimensional inhomogeneous Bethe-Salpeter Equation plays a key role for devising an equation for 
the Nakanishi weight function (a real function), as in the homogeneous case that corresponds to 
bound states and it has been already studied within different frameworks. In this paper, the whole 
formal development is illustrated in detail and applied to a system, composed by two massive 
scalars interacting through the exchange of a massive scalar. The explicit expression of the scat- 
tering integral equations are also obtained in ladder approximation, and, as simple applications of 
our formalism, some limiting cases, like the zero-energy limit and the Wick-Cutkosky model in the 
continuum, are presented. 

PACS numbers: 11.10.St,11.80.-m,03.65.Pm,03.65.Ge 
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I. INTRODUCTION 

The Bethe-Salpeter equation (BSE) (see also Q for an early and detailed review) is 
an important tool for investigating the non perturbative regime of many important issues 
in physics, within a relativistic framework. The huge tower of difficulties for obtaining solu- 
tions has been bravely faced with very different degrees of approximation and with different 
mathematical tools, e.g. by exploiting the classic Wick rotation trick jsj] and free propagators 
(see e.g. Ref. (4]), by three-dimensional (3D) reductions (see, e.g., Ref. js|) or by applying, 
in Euclidean space, the combination of Dyson- Schwinger and Bethe-Salpeter equations (see 
e.g. Ref. js] for a recent review). The peculiar field of application is the description of 
bound states, but one can find attempts for solving also the scattering-state problem (see, 
e.g., 0, 3Jo|). As mentioned above, the number of approaches for achieving solutions is 
quite large, and in this paper we focus on a particular scheme in Minkowski space. In order 
to avoid the well-known difficulties related to the cumbersome analytical behavior of the 
four- dimensional (4D) BSE, we adopt the framework proposed by Nakanishi in his seminal 
work [lo| , the so-called perturbation theory integral representation (PTIR) of i) the vertex 
function (that leads to the BSE for bound states), ii) the T-matrix (relevant for the present 
investigation of the scattering states) and iii) the general, multi-leg transition amplitudes. 
He carefully studied the parametric form of any Feynman diagram that contributes to the 
multi-leg transition amplitudes for interacting bosons, focusing on the denominator, since it 
governs the overall analytic behavior. 

Nakanishi was able, through a clever change of variable, to give a common form to 
any denominator appearing in the parametric expressions, and then he formally summed 
up the infinite series of Feynman diagrams. The final result, aside from the convergence 
problem of the perturbation series 10|, makes it possible to single out real, multi- variable 
weight functions (the Nakanishi amplitudes in what follows) modulating the denominator 
that generates the analytic behavior in Minkowski space of multi-leg transition amplitudes. 
It is noteworthy that a uniqueness theorem was proved for the PTIR transition amplitudes 



for bosonic systems 



io|. 



Summarizing, PTIR yields a well-defined framework where classes of approximations to 
the multi-le g tr ansition amplitude can be implemented in a systematic way, as shown, e.g., 
in Refs. ll|, ll2j , where the BSE was investigated first within a ladder approximation of the 
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kernel and then the cross- ladder class of diagrams was added. One could wonder if PTIR 
be able to deal with the analytic behavior of the non perturbative multi-leg amplitudes. 
For the scattering case, typical singularities like cuts, are certainly present in PTIR. For 
the bound states, one is interested in describing the vertex functions, which is related to 



the residue of the 
performed till now 



matrix at the bound-state poles. Moreover, numerical investigations 



11 



13| have shown that PTIR is able to accurately obtain both bound- 



state BS amplitude and corresponding masses. This strongly indicates [14J that solutions of 
the BSE, i. e. eigenfunctions obtained in a non perturbative framework, can be put in the 
form suggested by Nakanishi PTIR. In nuclear and hadronic physics, the Nakanishi PTIR 



has been applied to a wide range o 
state shown, e.g., in Refs. 



p 



rysical 



11-13 



problems, but so far restricted to the bound- 
15|-|20|. The relativistic scattering states play 



a relevant role in many hadronic processes, like in the light-meson decays of heavy hadron 
resonances (see, e.g., for a recent study [21]), or even in more exotic cases (see e.g. |22|). 
Therefore, the extension of the Nakanishi PTIR for solving the BSE in the continuum region, 
as one of the possible tools to deal with the 4D Minkowski structure without approximations, 
is appealing and deserves dedicated efforts. 

Recently, Carbonell and Karmanov proposed a new approach for obtaining, within PTIR, 



solutions of t 
fermionic 19 



re BSE with truncated kernel, for bound states of both bosonic [ill 1 121 



and 



20[ two-body systems, interacting through the exchange of a massive scalar 



boson (Yukawa Model). As in the case of Ref. 13[ (where only the bosonic case was investi- 
gated), they didn't make use of the Wick rotation, but, for the first time, they exploited the 



covariant formulation of the Light-front (LF) field theory 



obtaining a compact form for 



12j approximations. Within such 



the kernel of the BSE, in both ladder [ll| and cross-ladder 
an approach, it is possible to quickly get rid of the difficulties related to the analytic behavior 
of the BS amplitude in the Minkowski space, focusing the numerical efforts on determining 
a proper Nakanishi weight function, that remarkably are real functions, as mentioned above. 
The equation that determines such a weight function is directly related to the 3D integral 
equation for the valence component of the interacting-system state. In turn, such a 3D 
equation is obtained by projecting onto the null-plane the full BSE, after expressing the BS 
amplitude in terms of the Nakanishi PTIR. One may wonder how the valence wave function 



be able to generate the whole physical complexity of t 
a non perturbative framework, it formally turns out 



re full BS amplitude. Indeed, within 



24 



29j that by introducing the pro- 
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jection onto the null-plane, namely the analytic integration over the variable k~ = k° — k 3 , 
the full BS amplitude can be expressed in terms of the valence component and an operator, 
producing the richness of the Fock states on top of the lowest one. It must be emphasized 
that such a result has been obtained for both bound and scattering states of bosonic and 
fermionic systems. This gives us confidence that the PTIR approach, adopted for the bound 
state, could be extended to the scattering states. 

Aim of this work is to construct the integral equation that determines the Nakanishi 
weight function for the scattering states of two massive scalar bosons, exchanging a massive 



scalar (cf 



111 ] for the bound state). The main tool for obtaining the scattering equation is 



the projection onto the null-plane (LF-projection) of the BSE, as developed in a series of 



papers 



24j-[28|]. This method, not explicitly covariant, can be applied also to the bound-state 



case, obviously getting the same result like in Ref. [ll|, where the explicitly-covariant LF 
framework is considered. Moreover, the new integral equation is applied to the simple ladder 
approximation for the Yukawa model, and the extension of the Wick-Cutkosky model [si, [soj 
(a massless exchanged boson in ladder approximation) to the scattering region, is proposed. 

The paper is organized as follows. In Sect. HIl it is introduced both the general formalism 
of the BSE and the integral transform adopted for projecting both bound and scattering 



BSE's onto the LF hyper-plane (a short review on this topic can be found in Ref. [291] ) . 
In Sect. IIII[ a new integral equation for determining the Nakanishi weight function for 
4D scattering states is found, after briefly reviewing the bound state case. In Sect. IIV1 
explicit expressions for both the inhomogeneous term and the kernel in ladder approximation 
are given for the scalar Yukawa model with an exchanged massive boson. Within such a 
framework, some relevant cases are discussed in Sect. |Vj i) in subsect. IV A\ the zero energy 
scattering amplitude, ii) in subsect. IV Bl the Wick-Cutkosky model for scattering states 
and iii) in subsect. IV C\ a new form for the integral equation determining the Nakanishi 
amplitude for bound states, in ladder approximation. Sect. [VI] contains concluding remarks. 
In Appendices IA1 and IB1 more details on the formal developments of Sect. [IV] are given. 

II. GENERAL FORMALISM AND LF PROJECTION OF THE BSE 

The ingredients of our approach are in order: i) the PTIR of the BS amplitude, $(/c,p), 
and ii) the LF projection of the BSE introduced for obtaining a 3D equation for the valence 
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wave function, kj_), for both bound and scattering states. Once the 3D equation is 
obtained, one can determine the Nakanishi amplitude, that in turn it yields the 4D BS 



amplitude within the PTIR (lOj. It should be emphasized that the LF projection of the 
BSE in Minkowski space represents a smart way to deal with the singularities of the 4D 
once its analytic structure is known. 
The inhomogeneous BSE for a scattering state of two scalar particles with equal mass m, 
total momentum p = pi +P2, total mass M = \pp* and relative momentum k = (pi —p^jl 
reads as follows 

^ + \k^k l ) = {2^\k-k l ) + G ( i 2 \k,p) J ^iK{k,k\p)& + \k'iP,h), (1) 

where <3> < - + ) is the BS amplitude for the scattering state, % K, the interaction kernel, and fcj the 
relative incoming momentum. It is worth noting that the matrix elements of the T-matrix 
involved in the scattering process are the half-off- shell ones, since the incoming particles are 
on their mass shell. This implies that the square of the intrinsic four-momentum is given by 
kf = m 2 — M 2 /4 < (cf Sect. IIIip . The normalization of the free states is 

{k\q) = {2^{k-q) (2) 

and the free propagator of the two constituent particles is given by 

r (12)/r S _ r (l) r (2) _ i i ,o\ 

As well known, for a bound state one has the following homogeneous integral equation 
$ b (k,p) = GS 2 \k,p) J 0^iJC(k,k',p)Mk',p), (4) 

The above equations represent the starting point of our investigation. 

In the rest of the Section, the integral transform adopted for projecting both bound and 
scattering BSE's onto the LF hyper-plane is described. A reader already acquainted of both 
the Bethe-Salpeter formalism and the LF projection method can move directly to Section [TTTl 



A. The BS amplitude and the LF wave-functions 



As well known, the BS amplitude is defined as the matrix element between the vacuum 
(0| and a state \p), with total momentum p M , of the time ordered product of two Heisenberg 
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operators. For a two-boson system, one has 



$K,x 2 V) = (0\T{<PhW)<PhW)}\p). (5) 

In general, the state vector \p) can be taken in different representations, but within the 
LF quantization with only massive quanta, it can be written in terms of an infinite sum 



over Fock components (see, e.g. 31]), i.e. the Fock vacuum is an exact eigenstate of 
the full LF Hamiltonian (see also 32| for a recent discussion of the simplicity of the LF 
vacuum). The weights of the Fock states in the infinite sum are the LF wave functions 
(LFWFs), ipn/p, that have a probabilistic interpretation, namely the proper integral over 
|"0n/p| 2 yields the probability of the corresponding Fock state. It should be pointed out that 
the LFWFs provide an intrinsic representation of the composite system, as a consequence 
of the possibility to separate the global motion from the intrinsic one, since the LF boosts 
belong to the kinematical subgroup, and the mass operator is a Lorentz scalar. 

In terms of the four- momentum operator p = (p^ ,p~,p±), with jr 1 = p° ±p 3 , the invariant 



square mass operator reads (see e.g. |3l|) 



M 2 = jTp^ - p 2 ± (6) 
The eigenvalue problem to be solved is 

M 2 \p)=M 2 \p), (7) 

where Ai 2 is the eigenvalue. The eigenstate, \p), can be expanded in multiparticle Fock 
eigenstates {\ri}} of the free LF Hamiltonian, i.e a basis constructed by using the LF plane 
waves, The single particle free state describes the motion of a particle with mass m 
and it is an eigenstate of the LF momentum operators p + and p^, i.e. q = {<7 + ,q_i_}. The 
dependence upon the minus component, given by q~ = (m 2 + \<i±\ 2 )/q + , is understood in 
the notation we have adopted. The free-state normalization is 

(k\q) = 2k + {2n) 3 5(k + - q+) 5 {2 \k ± - qj = 2k + (2ixf 5 3 (k - q) (8) 

with the completeness given by 

d 3 q I a , | 
2q^f k> (9) 
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One should also recall that for a particle with mass m, one has 

d 4 qS(q 2 - m 2 ) 9(q°) = \ I d 3 q [ dq~ 5(q + q- - \q ± \ 2 - m 2 ) % H 



2 . 

% + ) ( 10 ) 

where the positivity of q + follows from the positivity of q°, since the delta imposes that 
q°= [g+ + (m 2 + |q ± | 2 )/g+]/2. 

Actually, each state \p) of the composite system is expanded in a complete Fock basis of 
non interacting n-particle states \n) with n > 2, since we are dealing with the BS amplitude 
of a composite two-boson system. Then one has 

oo „ 

\p) = 2 (2vr) 3 / M [ d ^] ^/ P (&,k iX ) |n; k^ + ^Ri), (11) 

n>2 

where the properly symmetrized free state with n particles are normalized, & = qf /p + and 
qi_L = kjj_ + £j pj_, with kjj_ the intrinsic transverse momentum of the % — th constituent, 
ip n /p(^i,'ki±) are the LFWFs (see below). It should be pointed out that {^,kjj_} is the set 
of intrinsic variables, fulfilling the following relations 

n n 

5> = 1 $> ± = (12) 

i=l i=l 

In Eq. (jlip . the phase-space factors are given by 



/[*]-n/*&(i-?>i. 



2(2*r)& 

/[^]-n/^(Ev). us) 

J 1=1 J v ; j=l 

Moreover, the kinematical nature of the LF boosts allows one to write the overlap 
(n\ £,iP + , kj_|_£jp_|_|p) in terms of the global motion and an intrinsic function, i.e. the 
corresponding LFWF, ip n / P (£,i, kjjj, viz 

n 

(n; £ iP + , k l± +Zip±\p) = 2p + (27rf5 3 (j2vi-p) Vty P (6,ki±) = 

i=i 

n n 

= 2(2vr) 3 5(l - 5^6) ^ (2) (E k ^) iM&ku) (14) 

i=l i=l 

For a bound system, the LFWFs are normalized according to 

£ J [<%i] [d 2 k l± ] \^j n/p (^k t± )\ 2 = 1. (15) 



Thus, given the Fock projection ( !T4|) . the state of a composite system can be determined 
by reaching through a suitable LF boost the intrinsic frame, where p_i_ = 0, and then 
solving the eigenvalue problem for the invariant mass operator M 2 , in order to obtain the 
set of LFWFs. It should pointed out that if one applies a dynamical LF rotation to the 
composite-system state, but still maintaining the Fock basis relative to the old orientation 
of the hyperplane, then the Fock content of the state changes. The content returns to be the 
same as before applying a dynamical rotation, if one adopts the Fock basis corresponding 
to the new orientation of the hyperplane. This is analogous to the Wigner transformations 
of the angular states, but generalized to the Fock space. 

Summarizing, from Eqs. an d ( II ip . one recognizes that all the Fock components with 
n > 2 are required for reconstructing the BS amplitude in Minkowski-space. This is easily 
understood, since the mass operator ([6]) is non-diagonal in the Fock space, due to the non 
commutativity of the interaction with respect to the particle-number operator. However, 
only the valence component (n = 2) is active in the calculation of the BS amplitude, (|5]), if 
we restrict the analysis to the equal LF-times case (xf = x^) (see below for details). 



B. The BS amplitude and its Fock content 



The general picture deduced from Eqs. ((5]) and 



take into account the results obtained in Refs. 



FTTj) . can be made more complete, if we 



24-: 



28 



for bound and scattering states of 



both two-boson and two-fermion systems. There, by using the LF projection technique (see 



also 



29|). it was shown that the valence component of the LF bound and scattering states 



contains enough information for reconstructing the corresponding BS amplitudes. It should 
be pointed out that a similar result for the bound state of both two-boson and two-fermion 
systems, has been obtained in Refs. [HI, 



12 



19], but by applying the Nakanishi PTIR to 



the BS amplitude. 

In what follows, the above mentioned relation between the valence component and the BS 
amplitude is briefly illustrated. From the field theoretical point of view, the states of the Fock 
basis, |n), are constructed by applying free-field creation operators to the vacuum state |0). 
Indeed, for a massive particle within the LF framework (let us recall that the quantization 
rules for the LF field theory 3]J are given at fixed LF time r = t + z/c 33[ ) , the eigenvalue 
of p + must be at least m, since the eigenvalues of p 2 and p° are both positive (the mass of 
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the particle yields the lower bound). Then, the vacuum, that has a vanishing eigenvalue of 
p + , cannot contain massive particles, given the conservation of the longitudinal momentum. 
This simple kinematical chain makes the Fock expansion meaningful (for a new vistas on the 



condensate in QCD see Ref. 32j). It should be pointed that the presence of the so-called 
zero-modes, that can be associated to particles with vanishing masses, spoils the triviality 
of the LF vacuum, i.e. p + \0) = does not more entail a vanishing number of particles. This 
means that the eigenvalue of p + is not sufficient for identifying the vacuum with respect to 
a degenerate state populated by zero modes. In this Section, we are considering a massive 
theory. 

In the composite system we are considering, the fundamental constituents appear in 
LF quantization as the excitations of the dynamical fields, <p(x), expanded in terms of 
bosonic creation and annihilation operators on the null-plane, with coordinates x~ = x° — x 3 
and x_i_, at fixed LF time x + = x° + x 3 31|. For x + = 0, the Heisenberg oper- 
ators turn into the Schrodinger one, and therefore the field can be written in terms 
of non interacting creation and annihilation operators as follows (notice that <-Ph{x) = 
exp{iP~x + /2) ip(x) exp(—iP~x + /2)) 

/rlh + rl 2 \< i / ~ \ 

where x = (x~ , xj_) and x ■ k = k + x~ /2 — xj_ • kj_ (in this notation x + = is dropped off). 
The creation and annihilation operators satisfy the following commutation relations 

a(k), a\k')] = 5(k + - k' + )5^(k ± - k' ± ) . (17) 



A one-particle free state is defined by \k) = (27t) 3 / 2 a/2/c + a^(k)\0), given the adopted 
normalization (cf Eq. OH])), and (x\k) = e . 

Let us express <&(xi,X2,p) in (jSJ) through its Fourier transform. Translational invariance 
imposes to $ the following form 

$(xi,x 2 ,p) = Hx,p) e~ ip - x , (18) 

where = (x± + x^)/2, x M = x± — x% and $(x,p) the reduced amplitude. It is related to 
its Fourier transform, <&(k,p), as follows 

d A k 



Hx,p)= I ™ e *~*(k,p) (19) 



where 

p"=p?+p£ F = ^^ (20) 

with p| 7^ m 2 (see Eq. (OQ)). The amplitude satisfies the homogeneous, (jlj), or the 

inhomogeneous, ([1]), BS equation depending if one is considering bound or scattering states. 
For X p = 0, one has 

®(x,p) = (0\T{cp H (x^2)cp H (-x>"/2)} \p) = 

= 6(x + ) (0\<p(x/2) e ~ iP ~ x+/2 <p(-x/2)\p) e ip ~ x+/4 + 

+ 6(-x + ) (0\ip{-x/2) e iP ~ x+ / 2 <p(x/2)\p) e~ ip ~ x+/i = 

= 9(x + ) ^e ip ~ x+ l\U\ip{x/2)\n') (n'\ e - iP ~ x+ /2 \n) (n\<p(-x/2)\p) + 

n,n' 

+ 6(-x + ) ^e- ip ^ +/4 (0|^(-x/2)|n / ) (n'\e iP ~ x+/2 \n) (n\tp(x/2)\p) (21) 

n,n' 

where \n) and \n') are states of the Fock basis used in Eq. ( fTTI) . It is easily realized that the 
matrix elements {n'\e~ tP x+ ^ 2 \n) are not diagonal in the Fock space, apart the case x + = 0, 
since the operator P~ contains the interaction in its full glory. Therefore, all the Fock 
components are acting in BS amplitude of the composite system. 



C. The BS amplitude and the valence wave-function 



The interesting; = leads straightforwardly 

component, ip n =2/ P , and the BS amplitude (see also 



to the relation between the valence 
23j|). In particular, one has 



lim $(x,p) = (0\ip(x/2) ip(-x/2)\p) 



x+->0+ 



(22) 



If we perform the 4D Fourier transform of the quantity <&(x,p) S(x + ) (i.e restricting the 
reduced amplitude to the LF hypersurface, or equal LF-times case (x + = xf — x% = 0)), 
one can extract the valence wave function ip n =2/p{^i k_i_). As a matter of fact, on one hand, 
one can write 

1 f dk~ f d 3 k 



d 3 x e-^ x <f>(x,x + = 0,p) 



2 J 2tt J (2tt 



$(k,p) / d 3 xe i{k ~^ 



dk 
~2^ 



(23) 
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On the other hand, one has 



d 3 x e %k ' x x 



2p\ 



d 3 x e~ il ~ x (0\ip{x/2)tp(-x/2)\p) 
d 3 q ± f d 3 q 2 



5 (P-Ql- Q2) ~J= ^n=2/p(£, kj 



2g+(27T)3 

2q+ J2q+ (2tt) 3 



2g+(2 



7T 



2p + (2vr) 3 x 



\a(Pi) a(p 2 ) a ] {qi) a j (q 2 )\0) 



p+y/2 



: p 2 + k) + (i-fc)- 



(24) 



Summarizing, the valence component is given by the Fourier transform of = + ,p), 

and it reads (cf 231 ]) 

^$(M = ^£(i-£)<M£,k ± ) (25) 



^n=2/ P (£, kj.) = ^= £ (1 - 



where ^ = £1 and we have introduced the notation 0lf(£, kj_) for future purpose (see the next 
subsection). The expression ( 1251) can be considered as a first bridge between the Fock state 
decomposition of a composite-system state, Eq. (llll) . and the BS amplitude. In particular, 
the integration over k~ projects the BS amplitude $>(k,p) onto the LF-hyperplane, i.e. 
x + = 0, and it leads to the valence wave-function. 

The LF projection can be directly applied to both the homogeneous and inhomogeneous 
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121 ] for the LF covariant approach for the bound state case). For instance, 



BSE (see Refs. 
from Eq. (jlj) one has for the bound state 



dk~ 



G<o 2) (k,p) [ ^ K(k,k',p)^ b (k',p), 



(26) 



2tt -"v--/ j 2vr ~ u x " ir ' J (2vr) 4 ' 

This mathematical step is a fundamental one in order to get an equivalent equation that 
allows one to determine the Nakanishi amplitude, eventually leading to the BS equation 
treated within PTIR framework lo| . (as discussed in detail in the next Section). In what 
follows, the LF projection is applied to both bound and scattering valence wave-functions. 



D. The LF projection of the BS equation 



The LF projection, applied to the BS amplitudes for bosonic and fermionic states, can 



be embedded in a more general framework 



24 J28 



where the relation (|25|) can be extended 



to the scattering BSE without relying on any perturbative approach: we reiterate that this 
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represents the main motivation of the present paper. It should be also mentioned that the 



general formalism of Refs. 



24-28 



allows one to fully reconstruct the 4D BS amplitude in 



Minkowski space from the 3D valence component (cf the following Eqs . (1391) and ( I40p ). 



25 



Let us briefly review the results for an interacting two-boson system 
for the generalization to a fermionic system and Ref. 



313,127 



23| (see Refs. 



for an introduction to the 



three-boson case). In order to get the 3D eigenfunction for the valence component from the 
4D BSE, the LF projection technique is combined with a Quasi-potential treatment of the 
BS interacting kernel (see 34|). This new ingredient makes it possible to single out from 
the 4D Green's function of the interacting system the "trivial" global propagation. Such a 
formal step is accomplished by using the LF projection, i.e. the proper integration over the 
k~ variable. Then, one remains with the valence dynamics described by a 3D effective LF 
interaction directly related to the BS kernel, through suitable projection operators. 

The main outcome for bound and scattering states is given by the 3D LF eigenequation 
for the valence component that can be obtained by applying the LF-projection technique to 
the BS equation. This can be rewritten as follows 



G-\p) \y)=[G \p)-V(p)] |¥) = 



(27) 



where an overall delta function for the total four-momentum conservation has been factorized 
out, Gq{p) is the free propagator and V(p) is the 4D interaction, obtained from some given 
Lagrangian (35)]. Their matrix elements are given by 



(pi\G (p)\p / 1 } = 5 4 (p' 1 - Pl ) G$*\k,p) 
(Pi\V(p)\p'i) = ilC(k,k',p) 



(28) 



where p = pi +p 2 — Pi + P21 k = (Pi — P2V2 and h! = (p[ — p' 2 )/2 . Then, following Refs. 



24H28| and the resume in Ref. 29j, one obtains the following 3D LF eigenequation (that 



holds for bound and scattering states, depending upon the boundary conditions) 



k±\g 1 (p) - w(p) \lp n =2/p) = 



(29) 



where the 3D free Green's function, go(p) describes the global propagation and w(p) is the 
effective 3D interaction, related to the 4D Quasi-potential (cf |34|). 



W(p) = V(p) + V(p)A (p)W(p) 
12 



(30) 



where the key ingredients is contained in A (p) = Gq(p) — Go(p) and it is represented by 
the auxiliary Green's function Gq(p), that is nothing else but the 4D image of the 3D go(p). 
Through Aq(p), one can factorize the intrinsic free motion, after subtracting the "trivial" 
global one. The 3D free Green's function is defined by 

9o(p) = ii G (p) lr (31) 

where the two new symbols 4- an d ir , that indicate a projection onto the LF hyperplane 
x + = 0, are given by 

el = J^- (K\ , U = / ^ \K) (32) 

Then, 3D quantities in Eq. ( 12 9 p have the following expressions 

go(py l =Mp) [GoipT 1 n (p) 

w(p) :=H (p) W{p) n (p) (33) 
where the /ree reverse LF projection operator and its LF-conjugated, are 

n (p) := Go(p) lr ^o(p)" 1 , Mp) := g^pY 1 ^G (p). (34) 

Finally, the 4D auxiliary Green's function reads 

G (p) = n (p) </ (p) n (p) (35) 

It has to point out that n (p) action is 3D — >• 4.D, while n Qo) acts in the reverse way The 
explicit expression of g^ip) is 

9o\p) = -i6(k+)6(p + -k+)kt (P + ~K) [p--Kn-Kn + *] ■ ( 36 ) 

where 

fe»l*l> = m2+ ,l kl± ' 2 l*l> (37) 

and the analogous for particle 2. The diagrammatic content (in the Fock space) of the 
Quasi-potential can be argued from the following formal solution of the integral equation 
that determines W(p) 

oo 

W(p) = V(p) + V(p)A (p)W(p) = V(p) ]T [A (p)V(p)} N , (38) 



N=0 
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where the Ao(p) insertions bring information on the intermediate particle propagations. 
A corresponding analysis can be performed for the LF effective interaction w(p), and it is 
physically quite transparent, since it can be carried out within a LF-time ordered framework. 
In particular, w(p) contains all possible LF-time ordered exchanges, corresponding to a sum 
over diagrams with any number of intermediate particles: this happens even for the ladder 
BSE. It should be pointed out that the ladder approximation of the LF eigenequation, Eq. 
([29]) . where w(p) is approximated by the first term in the power expansion in Eq. (138]) . does 
not account for the full complexity of the 4D ladder BSE. The physical reason lies in the fact 
that the iterations of the 4D ladder kernel and the 3D LF ladder kernel (time-ordered graphs) 
generate different intermediate states. The LF ladder kernel and its iterations contain only 
one exchanged particle (at a given global LF-time ) in the intermediate state, whereas the 
iterations of the ladder 4D kernel contain also many-body states, with increasing number of 



exchanged particles (stretched boxes, see e .g. [36(). This leads to a difference in the binding 



energies, which is however small 2J, |37|, |38[. In principle, for any 4D kernel given by a 
finite set of irreducible graphs, both BS ( 1271) and LF equation, Eq. ( 1291) . gives the same 
eigenvalue, once w(p) comes from the solution of ( 1381) . 

In the Quasi-potential framework the interacting LF reverse projection operator (see Ref. 



28|) 

n(p)= [i + A (p)w(p)} n (p) (39) 

leads to the following relations between the 3D valence component and the 4D BS amplitude 

|*> = n(p) \<t> LF ) (40) 

el |tt)= \<j> LF ). (41) 



The 3D valence wave- function, 4>lf, has been already introduced in Eq. (1251) . The analogous 
relations for the free case read 

|*o> = n (p) |^o) (42) 

4 |* > = |V>0>- (43) 

It turns out that the full complexity of the Fock-space in the 4D BS amplitude appears, 
not only through the effective interaction that determines the valence wave function, but 
also through the interacting reverse projection operator, H(p). 

14 



Finally, let us remind the 3D integral equations that follow from the eigenequation (129)) , 
for both the bound and the scattering states, viz 



(fc + ,kjso(pMp) \<f> b LF ) 

(k + ,k ± \l + g(p)w(p) |^o) 



(44) 



with g 1 (p) \vp ) = 0. 



III. INTEGRAL EQUATION FOR THE NAKANISHI WEIGHT FUNCTION FOR 
BOUND AND SCATTERING STATES 



In this Section, we illustrate how to obtain an integral equation for the Nakanishi weight 
functions for scattering states, starting from a BS equation for a system composed by two 
massive scalars, exchanging a scalar particle. The kernel is composed by the infinite sum of 



two-particle irreducible diagrams 



35 



and the self-energy contribution to the massive 
two-particle propagation are discarded, at the present stage. 

In order to have a suitable introduction to the scattering case, we first briefly discuss 
the S-wave bound-state case, within our LF approach. The same integral equation for the 
Nakanishi weight function has been devised by Carbonell and Karmanov 111, [12] , but within 



the explicitly covariant LF approach 



23]. 



The bound states 



In this subsection the S-wave bound-state integral equation for determining the Nakanishi 
weight function is presented by using the LF language introduced in the previous Section 



(namely a non explicit 



Nakanishi PTIR 
= i(-l) 



10 



y covariant formalism). The BS amplitude is written in terms of the 



ll| , gbipi' i z'\ k 2 ), as follows 
i Jo 

POO 

dz' I di 



1 



[7' + m 2 — jp 2 — k 2 — p ■ k z' — ie] 
g b {i, z'\ k 2 ) 



2+n 



[k 2 + p ■ kz' — 7' — k 2 + ie 



2+n 



where p 2 = M 2 is the invariant mass of the interacting system and 



m 



M 2 



(45) 



(46) 
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with m the constituent mass. For bound states k 2 > 0, while for the scattering ones k 2 < 0. 
n general, the power n in the denominator is a dummy integer parameter, and in Refs. 







111 . Il2j ] the minimal value n = 1 was chosen for the sake of simplicity. Bigger values of n 



13 



20] 



may result in a more smooth solution for 5^(7', z'\ k 2 ) 

The LF projection allows one to readily get rid of the singularity in Eq. (j45p . eventually 
obtaining the valence wave function, within the PTIR framework. In agreement with Ref. 



ll(, one can define ipb = V% ip n =2/ P as follows (cf Eq. flUD; 

dk~ 



£(i-0 



■Q b (k,p) = 
Zix 

g b {i, 1 - 2£; k 2 ) 



It should be pointed out that the denominator is always positive, and one can safely drop 
out the ie term. 

In order to strictly follow the notation of Ref. 
to the standard LF momentum as follows: 



(47) 



111 ], let us use the variables (7, z), related 



7 = *5 



1>z = 1-2£>-1 



(48) 



(1-z 7 



' , / 9b{j, z\ K 2 ) 

7 [7' + 7 + 2 2 m 2 + (l- z 2 )k 2 } 2 



(49) 



Then, one can rewrite 

4 Jo 

The integral equation for the Nakanishi weight function, ^(7, z; k 2 ), is obtained by inserting 
f )45|) in both sides of f l26|) . and using Eqs. f )47|) and ( l4T)]) (see Appendix A of [ill]) for details). 
One gets 

b(y, z\ k 2 ) 



H TT 



/o [7' + 7 + z 2 m 2 + (1 — z 2 )^ 2 ] 2 j 

where the kernel Vj,, appearing in the right-hand side of Eq. (150]) . is related to the kernel iJC 
in the BS equation, ([I]), as follows 



V b LF ( 7 ,z; 1 ',z , )= ip 



dk 



G ( 12) (k jP ) 



d 4 k' 



i)C(k, k',p) 



— ip~ 



dk~ 



(27r) 4 + p . fc/(Z , _ y _ K 2 + ze j : 

1 



d 4 k' 



2tt [(| + k) 2 - m 2 + ie] [(§ - k) 2 - m 2 + ie] 
ifC(k, k',p) 



x 



( 27r ) 4 [fc' 2 + p ■ k'z' -i -k 2 + ie] : 



(51) 
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Its explicit expression within both the ladder and the cross-ladder approximation for the 



kernel K, can be found in Refs. 



11 



12|. 



We emphasize that Eq. (I50p is equivalent to the initial BS equation ([I]), within the 
PTIR framework. Moreover, it should be pointed out that the mass, M, of the interacting 
system appears on both sides of Eq. (150 p and is contained in the parameter k 2 . Once the 
weight function (#,(7', z; k 2 ) is obtained, then we can reconstruct the BS amplitude (see Eq. 
( 145]) ). As already mentioned in the Introduction, the ability of the valence wave function 
to reconstruct the full BS amplitude is not only restricted to a perturbative framework, but 
it holds also in a non perturbative analysis, and even more for both bound and scattering 
]. This observation represented our motivation to face with the investigation 



states 



24 



28 



of the scattering sates. 

Finally so m e important remark, It turns out QQ that 9b{% ^) may be Z ero in 
an interval < 7 < 70, and that the exact value where it differs from zero is determined 
by Eq. ( l50l) itself. Moreover, the uniqueness of the Nakanishi weight function is ensured 

nn 



by a theorem (see ll0[). Following Refs. [l3|, Il4j . one can extract from ( )50l) the following 
eigenequation 

dr/ J ^ dz' V 6 ( 7 , z- 7', z'- K 2 )g b (i, z'- k 2 ) (52) 

For a discussion of the above equation within the ladder- approximation framework, see 
subsect. IVCl 



B. The scattering states 

The LF-projection technique, combined with the PTIR framework, has the appealing 
feature of systematically removing singularities, that appear in both BS amplitude and ker- 
nel, by integrating over k~ more simple analytical structures. We notice that the quadratic 
k° dependence in the propagators translates into a linear dependence upon k~ (see, e.g. 



40 



41|). Since within PTIR approach, the analytical form of any multi-leg amplitude is 
made explicit, the linearity of the pole structure attained through the LF technique is ad- 
vantageous. Then, it is natural to extend the investigation to the scattering case, increasing 
the complexity of the problem. 

According to the general treatment of a four-leg amplitude, as developed in Ref. 
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one can write the off-shell T-matrix as follows 

{k'»\T(p)\W) = U l=1J j\z t S (l- z )j j"^ x 

Q + (i,z t ) 



7 - [zi (p/2 + k'f + z 2 (p/2 - fc') + z 3 (p/2 + /c) 2 + z 4 (p/2 - fc) 2 + z b s + z 6 t + z 7 u] - te 

(53) 

where the Mandelstam variables are 

s = p 2 t = (k' -kf u = (k' + kf (54) 

and they satisfied the four-momentum conservation, that reads 

s + t + u = p 2 + 2(k 2 + k' 2 ) (55) 

Equation fl53l) has a redundant dependence upon the seven invariant squares, in order to have 
a compact form of the fully-off-shell four-leg amplitude (namely, without a decomposition 
in s, t and u channels). It is understood that the actual dependence is upon a set of six 
independent variables, given the constraint in Eq. fl55|) . 

Indeed, for evaluating the scattering states we need half-off-shell T-matrix, i.e. one has 
to put on their-own mass shell the two incoming particles with four-momenta p/2±ki. Then 
one has 



mT{p)\K) = n i=1 , 7 jT 1 dzi 5 y-~Yl **) /°° ^ x 



(56) 



7 - [21 (p/2 + fc) 2 + 2 2 {p/2 - k) 2 + 2; 3 m 2 + z A m 2 + z 5 p 2 + z e t + z 7 u] - ie 

where (p/2 + kj) 2 = (p/2 — kj) 2 = m 2 . As a consequence p ■ ki — 0, and therefore fcf = 
m 2 - p 2 /4 = k 2 < 0. 

Let us analyze in detail the ^-dependent part of the denominator in Eq. fl56l) . By 
exploiting the presence of the delta function in Eq. (!56j) for eliminating z 5 , it can be rewritten 
as follows 

z-y (p/2 + k) 2 + z 2 (p/2 - k) 2 + (z 3 + z 4 )m 2 + z 5 p 2 + z 6 t + z 7 n = 

p 2 

= (2- 2zi - 2z 2 - 2z 3 - 2z 4 - 3z e - 3z 7 ) — + (z x + z 2 + z 3 + z 4 + 2z e + 2z 7 )m 2 + 

(57) 



+ (^1 + z 2 + z 6 + z 7 ) 



(k 2 + V — - m 2 ) + C'p-k + (2k ■ k 
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where ( = (z 7 - z 6 )/(z 1 + z 2 + z 6 + z 7 ) and £' = (21 - z<i)j{z\ + z 2 + z e + z 7 ) belong to 
[—1,1]. Remind that the variables zi fulfill i) Yli=n z i = 1 an d ii) < Zj < 1. Furthermore 
by introducing a new variable 

,/ _ 7 + 2 ~ ?» 2 ) (zi + ^2 + -^3 + Zi + 2z 6 + 2z 7 ) _ p 2 (2-z 6 - z 7 ) 

[z\ + z 2 + z G + z 7 ) 2 (zi + Z 2 + Zq + z 7 ) 

that belongs to [—00, 00], recalling that (z\ + z 2 + z 6 + z 7 ) e [0, 1], then one can rewrite the 

half off-shell T-matrix as follows 

S + (7",C,C') 



(k"\T(p)\K) 



d( [ d£ i drf' 



7 



k 2 + v -j - m 2 + C'p ■ k + (2k ■ hi 



i58) 



— le 



where the weight function, Q + {i' 1C1O1 is t ne result of multifold integrations and changes 
of variables, as suggested by the previous calculation. The scattering wave can be obtained 
by considering the following vertex function 

i i 



(k»\G (p)T(p)\k?) 



1 rl PC 

d( / d£ / 

1 J-l J-c 



(I + k) - m 2 + re (§ - k) -m 2 + 



IE 



di' 



W,C,C) 



k 2 + £ - m 2 + C'P ■ k + (2k ■ hi 



— 7" + ie 



1 At f 1 ^ r a» (\ r ai * ^(t". ceo 

dQ I dQ J <ry / dct\ \ da 2 — 

1 J-i J-oo Jo Jo [D + ie] 



where the Feynman trick 
1 



pi— ai 

dcii / da 2 
Jo 



ABC Jo '""Jo '"" [a x A + a 2 B + (1 - a x - a 2 )C + ie] 
has been used. Hence, the denominator D is given by 

D 



(59) 



(60) 



1 - ai - a 2 ) 



k 2 + V — - m 2 + C'p ■ k + C2k ■ ki - 7" 



— h k ) — m 2 



+ a 2 



P uY 2 
k ) — m 



k 2 + — - m 2 + p • k z" + 2k ■ ^ z' - 7' 



(61) 



with z" = £'(1 — «i — 0J2) + «i — «2 and z' = Q(l — «i — 0:2) belonging to [—1, 1] and 
y = y(i — q/ 1 — a 2 ). Then, the vertex function becomes 

g^(i,z',z") 



(k»\G Q (p)T(p)\k?) 



dz' / dz" 
-1 J -1 



di 



k 2 + S- -m 2 + p- k z" + 2k- ki z' -'y' + ie 



(62) 
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where a factor (—i) has been inserted for convenience (cf H] for the bound state) and 



^+)( 7 ', z\ z") = i f da, f ^ da 1 Q + (-, -, - ) 
Jo Jo a \a a a 1 



with a = 1 — ccx — «2- The vanishing behavior of Q + for values of its arguments at infinity 
makes finite the twofold integration. 

A similar representation can be introduced for the scattering state b y using Eq. ( 1621 . 
with the same power as in the case of the bound state analyzed in Refs. Ill, 112], viz 

$(+\k,p) = (2n) A 5 {4 \k-k l ) + 



f 1 

I dz' / dz" 
'-i J-i 



d"f'- 



[7' + m 2 - \M 2 -k 2 -p-k z" -2k -hi z< - ie] ' 

{2Tif5^\k-k l )-i I dz' I dz" 

1 J -i j -00 

9^(7', z',z") 



d'i x 



(64) 



[i + 7 + k 2 - k~(k+ + fz" - fziz') - k+(fz" + krz') + 2z'cose^ru - ieT 
where we adopt the frame p = {M, 0}. Following the formalism of Ref. |l2|, Zi = 
—2kf/M (1 > \zi\ since the incoming particles have positive longitudinal momenta, i.e. 
p + /2 ± kf > 0), cosO = kj_ • kjj_, 7 = |kj_| 2 and 7 j = |kjj_| 2 . It should be pointed out that 
for the variable 7' a lower bound different from —00 could be possible, as discussed in Sect. 
ED (below Eq. (TTOOR 

The 3D LF scattering wave function is given by (tp^ = \[2 ip^} 2 , p ) 



ip^ +s> (z, 7, cosO) = p A 



'l-z 2 



dk~ 



$ (+) (&,p) 



/ 2) (27i) 3 S^Ck-k t ) + 



2tt 

4 

g(+\y,z',z;7i,Zi) 



dz' x 



[Y + 7 + 2 2 m 2 + (1 - z 2 )k 2 + f z z'(fzi + fcr) + 2z'cos6y/rTi ~ ^\ 
where the dependence upon the initial variables {7^ z^ in has been made explicit, for 
the sake of clarity. The result in Eq. fl65|) can be obtained by integrating over k~ the singular 
integral in Eq. (164]) as follows (cf Appendix A in Ref. [Ill ]) 
dk~ 1 



(65) 



2tt 



[Y + 7 + k 2 - k~(k+ + fz" - f Zi z') - k+(fz" + kr z ') + Zz'cosQ^rH - ie\ 



5(k + + fz" - fz iZ ') 



2 [y + 7 + k 2 - k+{f z" + fcry) + 2z'cos6y/YTi ~ ie ? 
i 5(—z + z" — ZiZ 1 ) 

M [7' + 7 + 2 2 m 2 + (1 - z 2 )k 2 + fz z'(fzi + k~) + 2z'cose^ffi - ie 



(66) 
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where z = —2k + /M, with 1 > \z\ (from the request of dealing with particles only, i.e. 
p + /2 ± k + > in the valence wave- function). 

It should be pointed out that the new term, k ■ k i: in the denominator of Eq. (1M|) does 
not produce substantial differences in the formal treatment that we applied to the bound- 
state case (cf Eq. (1451) ). Therefore, performing analogous formal steps, one obtains an 
inhomogeneous integral equation for the Nakanishi weight function, fr + ^(7, z, z'; 7i, z{). In 
particular, the inhomogeneous integral equation is obtained by inserting fl64l) in both sides 
of the following 4D equation, obtained from Eq. ([1]) 

[$«(M"(2^) 4 * (4) (*-**)] = G^ 2 \k,p) iK(k,ki,p) + 

+G^ 2 \k,p) J K(k,k\p) [& + \k',p) - (2ir)W\k> - h)] (67) 

and then by projecting onto the LF plane, as shown in Eq. (1261) for the bound states, viz 

f^-[^\k,p)-(2nr5^(k-h)] = f ^G$ 2 \k,p) i!C(k t k i ,p) + 
J 2n J Z7T 

: ~ G^ 2 \k,p) [ £L )C(k,k',p) - {2n)^\k' - h)] (68) 



2tt u v J (2tt) 4 ' 
After introducing Eq. (|64"|) in Eq. f )68|) and integrating over k~ , one obtains the following in- 
tegral equation for the Nakanishi amplitude, g^Xj, z', z; 7^, Zj), without angular momentum 
decomposition, viz 



1 [Y + 7 + z 2 m? + (1 - z 2 )k 2 + fz z'ifzi + k~) + Iz'coaB^fH - z'ep 

X LF (7, z\ ji, Zi, cos9) + 



+ / di d( dC'V^^ziyuZirfXX^^gt+WXX-nuZi)- (69) 

J-oo J-l J-l 

where the inhomogeneous term is given by 

/elk — 
——G { Q 2 \k,p) ifC(k,ki,p) = 
2tt 

= ~ P+ J ^[(E + ky-m 2 + ze] [(| - kf - m 2 + ie] %¥i ^ Kv) (70) 

and the kernel V S LF is related to the kernel iJC of the BS equation by (cf Eq. (1511) for the 
bound state) 

V s LF ( 1 ,z-- ll ,z l , 1 ',C,C,cos9)= ip + r^G^ 2 \k,p) x 



00 



2tt 



d 4 k" iJC(k,k",p) 

(2tt) 4 [fc«2 - + p ■ k" C + 2fc" ■ fci C - 7' + «e] 3 



(71) 
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It is worth noting that, modulo the value of k 2 , one formally has for the S-wave 

hm V s LF (-f, z; y u z h 7 , £ (', cos6) = V b LF (^, z; 7, C') 
C->o 



(72) 



Notice that the dependence upon 7$ and Zi is washed out by putting ( = 0. (see subset. IV CI 
for a discussion in ladder approximation). 



C. The scattering amplitude 

In order to complete the theoretical analysis, let us present the relation between the 
scattering amplitude and the Nakanishi weight function. For the sake of simplicity, the 
frame where the z — axis is perpendicular to the scattering plane has been chosen, then the 



scattering amplitude f(s, 9] 
matter of fact, one has (cf 



(with s = M 2 ), can be calculated from the BS amplitude. As a 



39|) 



f( 8 ,0) = T V(2^+ 1) e lSe sinS e P e (cos6 
lim (pi >3 / 2 | HC(p)\& +) ;p,ki) 



M 8n k'^k f 
1 



lim (tf,p\Go 1 (p)\& + hp,k i ) 



(73) 



re plane 
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28|), 



M 8lT k'->k f 

where k r = \Js/A — m 2 , k! = (p\ — p' 2 )/2 =, p\ + p' 2 = p and the orthogonality of t 
waves has been adopted (cf Eq. (OQ)). By using the LF projection method (Refs. 24, 
one can move from the 4D Minkowski space to the 3D LF hyperplane, simplifying the analytic 
integration, without approximations. Then, one can rewrite the scattering amplitude in 
terms of the LF 3D scattering wave function as follows 

1 



/M) 



M 8tt k'->k f 



lim {k'MG^\p)\^-p,h) 



— — Jim (k'\U (p) G^{p) Ii{p)\^pM 
M 8ir k'^k f 



where, from Eqs. ( HOI) and ( 142 jl . one has 

\ P ,k') = n \k') , ki) = n(p) \&, P ,~h) 



(74) 



(75) 



Recalling that (see Eqs. (133]) and ([39]), and Ref. |28|) 

n (p) Gu\ P ) n(p) = n (p) g^{ p ) [i + a (p) w( p ) } n (p) 
= n (p) Go 1 ^) n (p) = ^ 1 ( P ) 



(76) 
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since 



Mp) Go(p) -1 a (p) = n (p) - n (p) G , o(p)- 1 n (p)^(p)no(p) = o 



(77) 



one can write 



1 



M 87T fc'_>fc f 



(78) 



Finally, by inserting Eq. ( 136]) one gets 

/M) 



M 87r ( 7 ,z)^( 7/ ,2 / ) 4 



1 - z 2 



lim 



^ ( - + ' ) (z, 7, cos6 



(79) 



M 87r (7,«)^T-r/^/) V ~ 1 - z 2 

where ?/^ + ) is given by Eq. ( 165]) and the dependence upon 73 = |kjj_| 2 is understood. 

Equation (179]) can be written in terms of the Nakanishi representation by using i) the 
distorted part of tp^ and ii) Eq. fl69|) . viz 



/M) 



lim [7 + (1 - z 2 )k 2 + z 2 m 2 ] 



l LF h,z;7i,Zi) + 



(80) 



M 8 7T (7,z)->-(7/,a/) 

dy / rfC / dC' K iF (7, z; 7. *, 7, C, C, co^)^ + )(7', C, C; 7i, *) 

-00 7-1 7-1 

where 7/ = 7i and 2/ = ^. Notice that the factor 7 + (1 — z 2 )k 2 + z 2 m 2 vanishing for 
(7,3) — > (jftZf), is canceled out by the corresponding one in X LF and V S LF . This will be 
illustrated in the next section within the ladder approximation. 



IV. SCATTERING STATES IN LADDER APPROXIMATION 

In this Section, we present the ladder approximation of the integral equation ( 169]) . in order 
to determine the corresponding Nakanishi amplitude, g^fa', z', z; 7$, Zj). We also illustrate 
how uniqueness can be explicitly exploited in order to get a simpler integral equation, but 
with a more elaborated kernel. 
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A. The Nakanishi weight function integral equation 



In ladder approximation (starting from now, we drop the superscript LF to simplify the 
notation), Eq. ( IBTJj) reads 



1 fOO 



-i 7-oc [Y + 7 + z 2 m 2 + (1 - z 2 )k 2 + fz z'ifzi + + Iz'cosQjTfi - ie] 2 
= 7^ L \j,z;ji,Zi,cos9) + 

/oo pi pi 

H' \ d( d£V( L \ 1 ,z- nh z in "X,C',cose)g( + \i\tX'-,liiZi) ■ 
oo J-l J-l 

where is given by (see Appendix [X] for details) 

1 



(81) 



z]^ u z h cosV) = g 



[7 + (1 — z 2 )k 2 + z 2 m 2 — ie] 



,z;ji,Zi,cos9) (82) 



with 



e (L) (7 



9{z- Zi ){l-z) 



f3(z, zi) + 7(1 - Zi) - 2(1 - z)cos6> ^yfi - ie 
6(zi-z)(l + z) 



+ 



and 



P(—z, -z^ + 7(1 + Zi) - 2(1 + z)cos6 yfffi - ie 



M 2 



fi 2 + —(1 - z)(l + z^ - 2m 2 



P(z,Zi) = (l-z) 
In Eq. (IHTl) . (cf Appendix IB1 for details) is given by 

„2 

K (L) (7^;7i,^,7" ! C,C / ,cos 



+ (1 - Zj)m 2 



(83) 



(84) 



2(47r) 2 [7 + (1 - z 2 )k 2 + z 2 m 2 - ie] 



dv v 2 F(v,~/,z; -/",(,(', cost 



(85) 



where (dropping the dependence upon the external variables Zi, and 7, for the sake of 
simplicity) 



F(v, 7, z; 7 , C, C, cos #) = C(w, 7, 2;; 7 , C, C, cos#) + C(v, 7, -z; 7 , C, ~C, ~ z h cos 



(86) 



with 

C(v,y,z;-y",(,(',Zi,cost 



(1 + *) 



x 



0(C - ^ - *c) 



[ 7 + ^ m2 + K 2 (1 _ Z 2) + r>, z, ^ c, c, 7") + C, C; *) P?*** + zcose^/fn] - ie] ' 



(87) 
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In Eq. ( 1H7|) . one has 



X(v,z t ,(X') = v(l-v)(l + C-z t () 

r( Wi ,c,C',7") 



z(z,c,C;z l ) = C 



(1 + C'-*C) Ki-v) 

+ 



A/T 2 

C /2 ^ + ^ 2 (i + C 2 ) + 7" 



I) 



7 



(i + C'-*C) 



(88) 



Collecting the above results, one obtains the following integral equations for the Nakanishi 



amplitude in ladder approximation, g L 



(+) 



9l (i^' } z;j i} Zi) 



1 POO 

dz' / dj' 

i J-oo [7' + 7 + z 2 m 2 + (1 - z 2 )k 2 + z'zi + 2z'cos9y/yfi 



9 



h 

,;2 



- ttG (L) (j, ^; Ji, Zi, cosO) + 

1 — z z )k z + z^m z — %e\ 



+ 



2(4vr) 2 [7 + (1 - ^ 2 )k 2 + z 2 m 2 - ie] 



/oo pl rl 

H' \ d( / d£ x 



eft; W 2 7, z; 7", C, C, cosS) g£\i\ C, C; 7i, *i) (89) 

A solution of Eq. ( I89j) . obtained by retaining only the inhomogeneous term and assuming 
the uniqueness of the solution, is discussed in the following Section. Such a solution can 
yield some insights on the the analytic structure of the ladder approximation of the equation 
for determining the Nakanishi amplitude g^~\ i.e. Eq. (189]) . In particular, one could argue 
that a finite lower bound for the variable 7' could exist in order to reproduce the analytic 
structure of the inhomogeneous term. In particular, the inhomogeneous term is proportional 
to the global propagator 

1 1 



M 2 — 



4 ( m2 +^) _i_ ie " 7 + (1 - z 2 )k 2 + z 2 m 2 - ie 

- (1 — z 2 ) k 2 — z 2 m 2 . The lowest value is 7„ 



that generates a cut starting at 7 = — (1 — z 2 ) k 2 — z 2 m 2 . The lowest value is 7 m ,m = —m 2 < 
recalling that z 2 £ [0, 1]. Therefore, since 7 = |kj_| 2 , the lowest bound is zero. 

To conclude this subsection, it is shown the scattering amplitude, Eq. flHUj) . in ladder 
approximation,viz 



/ (L) (V 



■ 9 



M8vr 



g {L) (-f f ,z f ;-fi,Zi,cos9) + 



2% 2 



•1 i-i 
H' I d( / dC x 



1 J-i 



dv v 2 ?^, 7/, Zf\ 7", C, C'» cos °) 9 ( l\i", C, C; T*» z i) 



(90) 
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where Q^ L > and T are given in Eqs. ( 1551) and (155]) . respectively. 



B. Applying uniqueness to the integral equation 



Eq. ( 1891) can be rewritten in such a way that the denominator appearing in the lhs can 
be present in both the terms in the rhs. This allow to explicitly use the uniqueness of the 
solution of the integral equation and to obtain a simpler equation, as shown in what follows. 

By using the standard Feynman trick, one can rewrite X^ L > in a useful form for the 
following elaboration, namely (see Appendix [A] for a detailed discussion) 

Z i > L \ ll z ]ll ,z il cos9) = g 2 / dz' e(-z') / di x 



[7' + 7 + (1 — z 2 )k 2 + z 2 m 2 + ^-z z'Zi + Iz'cosQ ■Jiii ~ * e ] ' 
{6{z - Zi) 9[l-z + z\\ - Zi)] + 6{zi -z)6[l + z + z\\ + Zi )\} 



X 



(91) 



where 



la (z>) = z> (2k 2 - 



(92) 



As for Vs (L \ one can proceed through more subtle mathematical steps. This is thoroughly 
discussed in Appendix [B], here the final expression is given. One has 

,2 



V^ L \ ll z ]lh z hl "XX\cost 



di / dz'- 



g 



2(4tt) 5 



x 



J -1 [7 + z 2 m 2 + k 2 (1 - z 2 ) + i + z> (^zzi + 2cosd^ru) - ie] ' 
(1 + Z) 9(C-z-z t () Q'(z,z i - n ",i,z',C,C',^)+ 

6(z - (' + ZiQ Q'(-z, -z^ 7", 7, z', C, -(', /i 2 ) 



x 



(l + C'-O 

(1-*) 

(1 - C + ZiQ 



(93) 



where 



{z,Zi\l ,7 ,z ,C>C>^ ) = 6 



1 + c - *c 



„(zQ e(C)-fl(-*') g(-C) . 



A ( z, ,(,C;i",i';zi,iJ, 2 



(94) 
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where 



A (z, C, C; 7", V; *, /i 2 ) = £ J^(o) |^(o)) 
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%i(0)) (^(oM(c,c # ,r,* 2 ) -^MoMCc.c'.r,^) -/* 2 IJ (95) 

with y± (0) solutions of the following second order equation 

y 2 A{t, C, 7", k 2 ) + 2/ z u C, C', 7", y, /A A) + /1 2 = (96) 

The dependence upon z'/( > will become clear from what follows. 
In Eq. (I96p the coefficients are given by 

A/f 2 

^(C,C / ,7 // ,^) = C /2 ^ + ^ 2 (i + C 2 )+7 // 



4 

B{z, z , Zi, C, (, 7 , 7 , A* , A ) = A* +7 - 7 ^7 + A ^ - ^ (97) 



In Eq. f l94j) . one explicitly has 
1 



y±(o) 



2^(C,C / ,7 // ,^ 2 ) 

*<, C, C, 7", y , /i 2 , 0) ± s/&{z, 7, Zi , C, C, 7", Y, /i 2 , 0) - V ^(c, C, 7", k 2 ) 

= T (1 + C/ "f C) y±(0) (98) 

Notice that for A = 0, the dependence upon Zi in 2', 2^, (, 7", 7', /i 2 , 0) is dummy. 

Eq. (JM]1 can be simplified, noting that the first term, proportional to yf(0) 5(^(0)), does 
not contribute if fj, 2 7^ 0, while for /i 2 = (Wick-Cutkosky model, see the next Section) the 
second term, proportional to 0(yi(0)), is vanishing. 

By using Eq. fl9Tj) and Eqs fl93|) . fl89|) can be rewritten putting in evidence the analytic 
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behavior of all the terms, namely 



1 poo 

dz' / di- 



j^oo [i + 7 + z 2 m 2 + (1 - z 2 )k 2 + ^-z z'Zi + Iz'cosBjrfi - ie] 2 



g 2 / dz' e(-z') / di 



SW-z'(2k 2 -^ 2 )} 



[7' + 7 + (1 — z 2 )k 2 + z 2 m 2 + ^-z z'Zi + 2z'cos6 ^77^ — ie] ' 



{6>(z - 9[l-z+ z'(l - Zi)] + 9(zi -z)9[l + z + z'(l + Zi )}} + 



2(4vr) 



1 /"OO 



dz' / dj' 



d( j d(' I di' 

(1-z) 



[7 + z 2 m 2 + k 2 (1 - z 2 ) +-f' + z' {^zzi + 2cos6^ru) - ie]' 
+ 0(C'-s-*C) Q'(^;7",7V,C,CV)+ 

# ) (7 // ,C,C / ;7i,^) 



x 



L(i + c-ZiC) 
z - C + ZiC) 7", 7', ^, C> -C> K 2 ) 



(i-C + ^C) 

From the uniqueness of the solution of the integral equation (1991) . that we expect once the 
Nakanishi theorem for the vertex function (cf 10 1) is extended to the scattering case, one 
could write (see also [IJ ) 

g[ + \i, z', z- 7j , Zi) = g 2 9(-z') S[j' - z\2k 2 - /i 2 )] {9(z - z t ) B[l - z + z'(l - Zi )\ + 



(99) 



-e(zi-z) 9[l + z + z'(l + Zi)]} 



2(A-Kf 



1 rl 



I d( J d(' g { L\i',(,(';ii,Zi) x 



(l+c-ziO 

(l-z) 



9(C - z - ZiC) Q'(z, z % - 7", 7', z', C, C', V 2 ) + 
«<z - (' + ZiC) Q'(-z, -z^, 7", 7', J, C, -C, /i 2 ) 



(100) 



(l-C' + ^C) 

The first term in the rhs of Eq. (11001) yields the lowest order approximation to g^ ■ The 
following support for < 7' < fi 2 + 2\k 2 \ can be obtained by inspecting the delta function, 
combined with the constraint on z' . 



V. APPLICATIONS 

In this Section, some relevant limiting cases, i) the Nakanishi amplitude for scattering 
states at zero-energy (k 2 — > 0, fi 2 7^ 0), ii) the Nakanishi amplitude for the Wick-Cutkosky 
model in the continuum (/i 2 — > and k 2 < 0), and iii) a revisiting of the ladder kernel 
for bound states are illustrated in details. It is worth noting that at the end of the Wick- 



Cutkosky model subsection, it is presented a formal comparison between the kernel 



energy, obtained within our approach and the kernel that one can find in Ref. 
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or zero 



161 ]. for a 



S-wave bound state, putting a vanishing binding energy energy. It is rewarding to find a 
successful comparison with our formalism. Moreover, as shown in the last subsection, a new 
form of the integral equation that determines the Nakanishi amplitude for bound states, in 
ladder approximation, is obtained. This simple integral equation could suggest a different 
numerical investigation of the issue. 



A. The Nakanishi Integral Equation for the zero-energy scattering 



The choice z% = 7« = yields n 2 = — 7, — zfM 2 /4 = 0, namely a zero-energy scattering, 
given the physical meaning of k 2 = m 2 — M 2 / '4. For such values of Zi and 7$, one can simplify 
Eq. f )99|) . obtaining the Nakanishi integral equation for the weight function at zero energy. 
Notice that such a weight amplitude determines the scattering length through Eq. ( 1901) . 
After inserting z% = 7$ = k 2 = in Eq. ( 1991) . one gets 



1 dz' r d ^L\j^^i=zi=G) 

-1 J-00 h' + 7 + z 2 m 2 - ie] 2 

1 poo 

9 2 



dz' e(-z') / di — 



5{i + z'^ 2 ) 



j-00 [7' + 7 + z 2 m 2 - ie} 2 
{9{z) 9(1 - z + z') + 9{-z) 9(1 + z + z')} + 







2(4?r) 2 
d( [ d£ 



dz' 

1 

00 



00 



[7' + 7 + z 2 m 2 H — ie]' 



x 



cZ 7 " 



:i + z 



i + C 



jr9{C-z) Q'(z,Z i = 0;i',i,z',C,C,^)+ 



£ 0(z - CO Q'(-^, * = 0; 7", 7, C, -C, /i 2 ) 



^ +) (7 // ,C,C / ;7, = ^ = 00101) 



One immediately realizes that the integration over z' can be performed in the lhs, introducing 
the suitable S-wave Nakanishi amplitude for the present case. Then, one gets 

9^(i,z) 



00 



H-rj 



[7' + 7 + z 2 m 2 — ie} 2 
9(i) 



9- 



d^f' 



A 4 J -ex 

9 



9 



2(4vr) 



[7' + 7 + z 2 m 2 — ie\ 
dj' 



[0{z) 9(1- z- i/^ 2 ) + 9(-z) 9(1 + z- 77/i 2 )} + 



[7' + 7 + z 2 m 2 H — ie\ 



d(' / di' 



'1 + z 



'jrOtf-z) 7W,Y,CV) + 



;i+c 



(102) 
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where 



9ol\l', z) = I dz' g^'il', z', z; 7i = z { = 0) 



dz' 5{i + z'fi 2 ) = — 6{i) 



i 



^ ) (7",0= / dC 9 r'(i"X,C; 7i = ^ = o) 



(103) 



-i 



and the possibility to integrate over ( the function , (7",C)C';7t = Zj = 0) is related to 
the independence upon £ of the function 



<-7-/ // / >/ 2\ 

As a matter of fact, from Eq. ()94p one has 

T(z, 7 ",7',cV) = 



y dz 7 Q'(z,Zi = 0; 7", 7', z', C, /^ 2 ) • 



1 dz'fl + z 




1 



^,C',7',7'V 



z' Vi + c c 



[0(z') (9(C) - 0(-z') e(-C)] x 



where 



(l±f) z 
(i+C 7 ) 



9 



(104) 



5(^(0)^(0) -%,(0)) 



2/i 2 j/i(0) 



x 



y?(o)^o(C',y)-/* a J 



(105) 



with the dependence upon z'/( (always positive !) explained in what follows. From Eqs. 
(1971) and (1981) . one can explicitly write 



M 2 

A(C,7) = C —r + l =C "2 +7 



z', 0, c, c, y, 7", ^ o) = /i 2 + 7" - y| 

y±(o) 



2A(C',7") 



(V + 7" - y£) ± y (V + 7" - y^) - va(c, 7") 

(i + C 1 ) y±(o) 



d . s 

TTT-^(O) = =F n 



(106) 



(/x 2 + y-7'^) 2 - VA(C',7") 
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This allows one to single out the dependence upon ^ of the integrand in Eq. (11041) and, 
after introducing the variable 



z' (1 + 
X = — 



one has 



C (1 + *) 

1 dx 



x 



An Z,X- 



+ c v y u 2 



(107) 



This completes the proof that T does not depend upon (. Finally, from uniqueness, one can 
rewrite the inhomogeneous integral equation ( 11021) in a simpler form, viz 

z) = ~ 2 0{i) [*(*) 6(1- z- Hf) + 6{-z) 6(1 + z- i/f)] + 



2(4vr) 2 
'l-z 



d(' / di' 



'l + z 



i + C 



£*(*-C) T(-z,T>W,-C,rt 



(108) 



(i-C) 

Such an equation could have a direct application in the study of the relativistic effects of 
the scattering length. 



B. The Nakanishi amplitude for the Wick-Cutkosky model in the continuum 

The Wick-Cutkosky model 0, [k]], namely two massive scalars interacting through a 
massless one (/i 2 = 0) in ladder approximation, can be extended to the scattering case. The 



integral equation for the Nakanishi amp 
far only for the bound state case (cf [2j, 
becomes 



itude, for t 
3j, ydl and 



lis widely adopted model, was known so 



161 ] for a LF approach). Equation (jlOOp 



(V, z', z- 7i , Zi) = g 2 6(-z') 5{i - Iz'k 2 ) [6(z - Zi ) 6[1 - z + z'(l - Zi )\ + 
+6( Zi -z)9[l + z + z'(l + z^]] + 

2 roo rl rl 

dl" / d( / dC # ^(7*,C,C # ;7<,2i)x 



2(4vr) 



) 2 

[l + z. 



(i + C'-^C) 
(l-C' + ^C) 



6(C - z - ZiC) Q'(z, Zl - 7 ", 7 ', C, O /i 2 = 0) + 

6( Z - c + Z(C) Q'(-2, 7", 7', ^, C, -O /i 2 = o) 



(109) 
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where 

V 1 + C - ZiC J z> 

£*(0) %o(0)) (110) 
with |/o(0) being the non trivial the solution of the following second order equation 

A((, C, i")y 2 + yB(z, z', zt, c, (', y , /i 2 = o, A) = o (in) 

Notice that the trivial solution, y = for any A, does not contribute, due to the presence of 
dy (0)/d\. In Eq. f illip the coefficients are given by 

^(C,C / ,7") = C' 2 ^ + ^ 2 (i + C 2 ) + 7 // 

B(z, z\ C, C', 7", ^ = 0, A) = 7 " - yl + A (1 + C ^~ f C) (112) 

2/ (1 + Zj 



Then, one explicitly has 



Z/o(0) 



^(c,c,y',« 2 ) 

It has to be pointed out that, if A(C, C'? l") — 0, one has only the trivial solution y = for 
any A. 

Equation ( I109[) can be rewritten as follows 

2i+k7' , z', *; 7i, = <? 2 O(-z') 5{i - 2z'k 2 ) [6(z - zt) 6[l-z + z'(l - Zi )\ + 
+9(zi-z) 6[l + z + z'(l + z l )}} + 

6»(C - * - *C) f r \+ z , - t)+ 0{z -(' + ZiC) 1 



92 — r ^7^r<#fe4^C,C';.,,)x 



2(4ir) 2 * ' 7-oo i-i i-i .4(C,C',?",k 2 

- 2 - *> 6 (iTc^ic - ?) + e(z - c + z '° 9 (l^cfic - f ) H 

Notice that the extrema of the integration on ( in the second and third terms are different. 
It is easily seen that, in the rhs, the positive z' decouples from the negative one, obtaining 
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a homogeneous integral equation for O(z') g^{l' ■, z', z; ji, Zi) and an inhomogeneous one for 
O(-z') g^r{^' ', z' ', z;^i, Zi). Therefore, the Nakanishi amplitude for scattering states within 
the Wick-Cutkosky model is given by the following inhomogeneous integral equation (it is 
inhomogeneous for matching the boundary condition for scattering states) 

gtwii, z', z; 7i , zt) = g 2 0(-z') 8{i - 2z' k 2 ) [9(z - z % ) 9[1 - z + z'(l - z t )] + 
+6(z t -z) 9[1 + z + z'(l + Zi)]] + 



9' 



2(4tt) 2 v 7 i-oo 7-i J-i A(t,tw,K 

- « - -0 - + - C + ,0 * - 9 (H5) 

Let us assume the following separable form for 

(Y, z', *l 7<, zj) = ^(-^) - 2*'« 2 ) 4U(^'> *l 7i, • 
Then, by integrating both sides of Eq. (I115P on 7', one has (recalling that (/z' > 0) 

h^(z', z- 7i , Zi ) = g 2 {6(z - 9[1 - z + ^(l - *)] + 6{ Zi - z) 9[1 + z + z'(l + *)]} + 
ff 2 /"° f 1 , t , 1 7 (+) / >- >-/ x 

2(4*y L< c*Mi + K m + Q2 h w ( ^^" H}y - 
{ - * - *> * (t^c - ?) + 9(z - c ' + 2 '° 9 (to? " ?) } (116) 

It is worth noting that the solution h^y is perfectly compatible with the separable form we 
assumed, and that C' 2 ^f~ + k 2 (1 + C) 2 = 2(k 2 , k 2 ) is different from zero, as discussed 

below Eq. {USD- 
Summarizing, the solution of Eq. (I116p . together with the separable form mentioned 
above, allows one to obtain the Nakanishi amplitude for scattering states within the Wick- 
Cutkosky model in the continuum. 

As a formal check of our previous elaboration, in what follows we devise the same kernel, 
one has for the Wick-Cutkosky model for bound states in the limit n 2 — > 0. This allows us 
to make a successful comparison with the result for the S-wave, in the limit of zero-binding 



energy, that one can get from Ref. 
In the limit 



16] 



2 2 ~ ' z , r\ 

K = 771 — = —7,- — Z: — > 



M 2 „ikP 
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Eq. (I116P becomes (dropping the dependence upon z% = ji = 0) 

h^(z', z) = g 2 [9(z) 9[l-z + z'\ + 6{-z) 9[1 + z + z'\\ + 

{ e(C - z) e\z> - c^±^\ + «( z -OV-C^|]} ("7) 

By integrating over z' (let us recall that > z' > — 1), in order to match the S-wave 
dependence, one can write 

= PC*) (1 - *) + (1 + *)] ~ ^y 2 J\ f £ dC ^ OC, CO x 
° dz' ( 6»(C' - *) - C + 0(z - 9[z' - C ' 



(i + C) J v * (i-o J 

g 2 [9(z) (1-z) +9(-z) (! + *)] + _£_ x 

>^ { w _,)fi±i + *(.-og^} -Sko (118) 

where < (1 ± z)/(l ± CO < 1 and 

4U(co = /Voc,co (H9) 

The S-wave Nakanishi amplitude is even for z — >■ — z transformation. 

It should be pointed out that the inhomogeneous term has the same form, when the 
zero-energy limit is performed first, and the limit fi 2 — » 0, as briefly shown below. As a 
matter of fact, from the inhomogeneous term in Eq. (1108j) (corresponding to n 2 — > 0), one 

gets 

,2 



9 

,2 



hm ^ (9(1 - z - tV) + *(-*) «(1 + z - i'/fu 2 )] 



= ^(y)[^)(i-z) + e(-«)(i + z)] (120) 

where the the following representation of the delta function has been used 

hm -V(7) [//(l - \z\) - i] = (1 - \z\) 5( 7 (121) 
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Such a relation can be deduced by i) introducing the infinitesimal quantity (1 — |z|)/i 2 , ii) 
folding the previous limit with a generic function /(Y) and iii) integrating over 7'. 

As to the kernel, as above anticipated, it coincides with the one given in Ref. 16| . 
for a bound S-wave state, but in the limit of zero binding-energy. This comparison, even 
if it does not apply to a physically meaningful case, due to the presence of an essential 
singularity, ( >2 m 2 , gives us confidence in our formalism, since it is able to reproduce the 
result corresponding to a well known case, after performing the non trivial limits fi 2 — > 
and then k 2 — > 0. 

Finally, it is important to remark that even the non relativistic Coulomb problem, for pos- 
itive energies, studied through the Lippman-Schwinger equation, has an infrared divergence 
that must be treated carefully. In the Wick-Cutkosky model this problem is present and 
has to be addressed to provide meaningful solutions for the relativistic case in the scattering 
region. 



C. The S-wave bound state in ladder approximation revisited 



The integral equation that determines the Nakanishi amplitude for a S-wave bound state, 
in ladder approximation, can be rewritten in a simplified form, with respect to the one 



presented in Ref. 



ll| , once the uniqueness of the solutions is exploited. In particular, one 



can start with the expression of Vs L ^ given in Eq. (I93il . carefully performing the limit ( — > 
and taking a negative value for k 2 in the final result. 



From Ref. 



Ill ], one has for a S-wave bound state 

9b(l',z;K 2 ) 



00 

[7+7 
00 pi 



z 2 m 2 + (1 — z 2 )k 2 



ie\ 



di> / <H(7,^C'W,CV) 



122) 



A direct comparison between Eq. (IBlj) and Eq. (IB2j) . modulo the value of k 2 , leads to the 
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following relation for the interaction kernel VS 



V b iL) (7, z] l\ CO = Jim V S (L) (7, z; 7^ z h 7", C, (', cost 



2(4vr) 



lim 



di I dz' ; 

-1 [7 + z 2 m 2 + k 2 (\ - z 2 ) +i + z' {^z Zl + 2cos9^TU) - ie] ' 



X 



(1+*) 



(i + C'-*0 

(1-z) 



9(C - z - Zl Q Q'(z, Zl - 7", 7', z', C, C, ^)+ 

0(* - c' + ZiC) Q'(-z, -Zii 7", 7', C, -C, /i 2 ) 



(123) 



(i-C' + ^C) 

where <2' can be found in Eq. ( 1B37I) . 

Indeed the limit on ( must be carefully performed, since Q' contains two terms, in the first 
one there is 8(() and in the second one there is 6(—Q. Therefore the limit in the two contri- 
butions must be ( — > + and C — > CL, respectively. In what follows this is understood. The 
relevant limit to be investigated is (the one corresponding to Q'(— z, —zf, 7", 7', z', £, /1 2 ) 
can be carried out analogously) 



nil a 1 /-i 2 2\ 



dz'- 



(1 + *) 



lim 

Q'(z,z z]1 ",j',z',(,C,V 2 ) 



eiC-z-ZiC) x 



[7 + z 2 m 2 + k 2 (1 - z 2 ) +j' + z' {^zzi + 2cos6^/Tfl) - ie] 

lim 7 - 0(C' - ^ - ZiC) x 

C-K) (l + C'-^C) 

1 A (z, C, C, 7", Y, 7i, 2 ) 

1 ^ [7 + 2 2 m 2 + k 2 (1 - z 2 ) + 7" + z> (^j-zzi + 2cos#^77^) - ie] ' 
l + z 



x 



c-o i + C'-*C 



- t W) <KC) - o(-z') 9(-c)\ 



e(C - z - Zi Q) e(c) 



c fix 



1 + C - ztC 



X X 



[7 + z 2 m 2 + k 2 (1 - z 2 ) +Y + x( (^zzi + 2cos9y/fyi) ~ ie]' 



lim, < 1 + «> 

C^O (1 + C - *C) 



X X 



A(z,x,C,C / ,7 // ,7 / ,^,7i,/^ 2 ) 



[7 + z 2 m 2 + k 2 (1 - z 2 ) + 7' + xC (^z* + 2cos0^/77~) - ie] 
where the change of variable x = z'/( has been inserted. By recalling that 

1 l + z 

> 



id - l+C'-^C 



(124) 



(125) 
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since 1 > |£| and 9{Q' — Zi( — z) one can write 

(1 + *) 



•/() 



iC da; 



C^o (1 + C - ZiC) 



2 (126) 



x [7 + z 2 m 2 + k 2 (1 — z 2 ) + 7' + [^-zZi + 2cos6^yjjij — ie\ 
where the following identity has been exploited 



"° dx 

:-o / - 



1 X 

c 



,'ici dx 
<) - 



X 



(127) 



Then, one gets 



(l+f) fl(C - g) 

(1 + CO [7 + 2 2 m 2 + K 2 (l-2 2 )+y-ze] 



X 



^^A(^,C' )7 ", 7 > 2 ,, 2 ) 
x 



(128) 



where 



A (2, x, (', 7", 7 , /i 2 , k 2 ) = 



i=± 



x 



yf(0) MC, 7", « 2 ) - A* 2 



(129) 



with 



y±(o) 



X 



2A(C',7",^; 

B 6 (*, x, C', 7", V, aA 0) ± ^ B 2 (z, z', Zi , C, C, 7", 7', /i 2 , 0) - 4A(C, 7", ^ 2 ) 
0* (1 + *) C, 7", Y , A* 2 , 0) - VA(C> 7", ^ 2 ) 



(130) 



and 



A(C / ,7 // ^ 2 ) = C /2 ^ + ^ + 7" 

10 / J-l II I 1 r\\ 2 1 // 7 

06(z,£,C,7 '7 ^ ,0) = /x +7 

X 



M 2 



(131) 



The positivity of 7" in Eq. (1122ft leads to the positivity of At- This implies that £?& be 
negative (cf the end of Appendix IB"]) , and eventually 7' be positive, consistently with Eq. 
(11221) . analyzed in Ref. (ll|. This could also be expected on physical grounds, since the 
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particle production, or a cut in the mathematical language, must be avoided in the bound 
state kernel V b (7, z; 7", CO? namely in the denominator appearing in Eq. ( 1 1 3 2 p . Therefore 
the interaction kernel becomes 



J2 roo 1 
^)( 7 ,,;/,0 = --|— / H ; X 

lykKy J [7 + z 2 m 2 + k 2 [1 — z 2 ) + 7' — zej 



jt— 7a 0(C - 2)/ — A (z,x,C,7 ,7 ,A* .« ) + 
I 1 + s J Jo x 



+-^r|y ^ - CO 1^ ^ A (-*, x, -co 7", Y, a* 2 , * 2 ) } (132) 



The new form of the integral integration for 5^(7', z; k 2 ) is given by 

2(4^ 2 



2 /"OO I>1 

V-1 



X 



'1 + 2;) „ , u . 1 fT+C 7 r 



Tf--^ o(c -z) - a (z, x, co y , y , ^ ^ 2 ) + 

(1 + C) Jo ^ 

©(^ - co jf ^ f a (-,, x, -co y , y, ^ k 2 ) | (133) 



VI. CONCLUSION 



We analyzed the Bethe-Salpeter equation for scattering states by exploiting the Nakanishi 



perturbation theory integral representation of the multi-leg transition amplitudes [10|. In this 
way, one has the possibility to single out their analytic behavior in Minkowski space and to 
exactly project (namely explicitly integrating over the variable k~) any multi-leg amplitude 
onto the null plane, so that the relative Light-front time (between free legs) can be eliminated. 
In this context, we have addressed the BSE for the scattering states, in Minkowski space, 
extending the work of Ref. [lj], [ijj], where the explicitly-covariant Light-front approach 
was applied to obtain bound states for a massive two-scalar system, interacting through a 
massive scalar exchange. 

A key ingredient of our work is the one-to-one correspondence between the BS amplitude 
and the LF valence wave function of the interacting system. Such a relation is implemented 
through an operator that is able to produce the full complexity of the Fock space on top of 



the valence component 24j-|29l| . and noteworthily it makes feasible the study of the Nakanishi 



amplitude through the valence component, without loosing any physical content. 

38 



The LF projection technique allows us to write down an exact 3D equation for the valence 
component from the 4D BSE, both for bound and scattering states. In turn, given the re- 
lation between BS amplitude and valence component on one side, and the Nakanishi weight 
function on the other side, one can determine the last one, through a more simple math- 
ematical treatment, namely by using the 3D equation for the valence component without 



facing with the complexity of the 4D Minkowski space (see for the bound case [13] ). 

In the case of scattering states, it has been necessary to consider the PTIR of the half- 
off-shell T-matrix, i.e. a four-leg amplitude (while, for the bound state, it is needed the 



PTIR of the vertex function 



11 



12j|). Then, applying the projection method, the four- 
dimensional inhomogeneous BSE has been exactly reduced to a three-dimensional equation 
for the valence component, that in turn it has allowed to devise an equation for the Nakanishi 
weight function (a real function), without angular momentum decomposition. The relevant 
equation, for a massive two-scalar system, interacting through a massive scalar exchange in 
the continuum, is Eq. ( |69|) . that for the sake of clarity we report in this Conclusion, viz 

\,r H s <+>(y,^;7„*,) 



i J-oo [Y + 7 + z 2 m? + (1 - z 2 )k 2 + fz z'(fZi + + 2z'cose y /Tf i - ief 
X iF (7, z; j h z h cosO) + 

oo pi pi 

di d( dC'V s LF ( 1 ,z; 7l ,z i , 1 ',CC,cose)g^\ 1 'XX , ;i l ,z l ) 



where the inhomogeneous term, X LF , and the kernel, V S LF , are given in Eqs. (!70|) and (!7T|) . 
respectively. As a by-product, we have obtained the scattering amplitude in terms of the 
g^ + \ X LF and V S LF . This relation will be useful for phenomenological studies. 

The explicit expression of the previous integral equation has been also obtained for the 
4D BS kernel in ladder approximation, preparing the matter for forthcoming numerical 
investigations. In particular, a simpler integral equation for determining the Nakanishi 
amplitude, in ladder approximation, has been worked out by explicitly applying uniqueness, 
see Eq. ( 11 00p . It is worth noting that, in ladder approximation, the validity of the uniqueness 
in the continuum can be checked by obtaining the Nakanishi weight function, through the 
two possibilities given by the integral equation in ( 1991) and the one in (llOOp . Moreover, it 
should be emphasized, from one side, the simplicity and the benefit of the LF method for 
obtaining the kernel of the integral equation for scattering states, and from the other side 
the elaborated mathematical steps necessary to explicitly use uniqueness in order to simplify 



the form of integral equation (1991 (see e.g. the work of ref. jl3| ) 
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We explored our formalism in ladder approximation, by investigating two limits: i) the 
zero energy limit n 2 — > and ii) the Wick-Cutkosky model in the continuum, that corre- 
sponds to an exchanged massless scalar, namely fi 2 — > 0, obtaining for this case a separable 
form for the Nakanishi weight function. The cross-ladder contribution will be considered 
elsewhere. Moreover, a new form of the ladder kernel for bound states has been provided, 
allowing a simplified form for the integral equation determining the Nakanishi weight. 

Nakanishi conjectured 1^] that his approach could be valid beyond the perturbative 
regime. Indeed, numerical solutions of the bound-state problem have shown that, for trun- 
cated kernels, his method can be successfully applied. The investigation of the validity of 
the Nakanishi approach to non perturbative problems in the continuum, as a reasonable 
extension of what has been already done for the bound-state poles of the transition matrix, 
can be performed by solving the integral equations obtained in this work. 

Finally, the Nakanishi PTIR, applied to bound and scattering states, opens the possi- 
bilities of studying many new issues, since it is not constrained to 3+1 dimensions or to 
two interacting particles. For instance, let us mention that it could be applied to three- 



interacting bosons (see e.g. 



211 ] for the three-boson BSE and the investigation of the final 



state interaction in three-body decays of heavy mesons), or to consider the extension to 2+1 
dimensions and fermionic systems (useful, e.g., to treat the Dirac electrons in graphene, for 



reviews see 



12 



43|). 
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Appendix A: The inhomogeneous term in ladder approximation 



The inhomogeneous term, present in the 3D LF equation for scattering states (cf Eq. 
(J55]l ). is given by 

X iF (7, z;ii,Zi,cosO) = p + J ^-Gq 2) (k,p) ilC(k, k h p) = 

= fp+ I tr [(| + kf -mH ie] [(§ - k) 2 -mH ie] ^ Kp) (A1) 
Let us calculate the term in ladder approximation. Within such an approximation, the 
kernel z/C is given by (cf Ref. 

where fi is the mass of the exchanged scalar (g 2 = 16irm 2 a) and the momentum transfer 
is given by p\ — pu = (p/2) + k — (p/2) — ki. Let us recall that ki = (pu — P2i)/2 and 

P = Pli+P2i = Pi +P2- 

Then, in a reference frame where pj_ = and p^ = M, one has 

X (L) (7, z\ 74, z h cosQ) = ig 2 p + x 

dk~ 1 1 1 



2tt [(£ + k) 2 ~m 2 + ie\ [(§ - k) 2 - m 2 + ie] (k - h) 2 - /1 

. 2 + f dk~ 1 

W P J ^ (M/2 + fc+) (M/2 - fc+) (jfc+ - jfef) X 

1 1 



2 -l ie 



[(§ + fc)~ - (1 + fc)~ + te/(M/2 + fc+)] [(1 -*)--(§- + ie/(M/2 - fc+)] 
1 



x 



dk 



2tt AP(1 - z 2 ) (z - Zi) 

1 1 



X 



f + k- - (I + fc)~ + <2e/[M/(l - *)] f - fe" - (1 - fc)~ + <2e/[M(l + *)] 

fc- - fcr - (k - h) m - i2e/{M(z - Zi )] (A3) 
with z = -2k + /M, Zi = -2kf /M = 2k~ /M (since p ■ ki = 0) and 

f£ + frY- = - 2(m 2 + 7) 
l 2 jon Plon M(1-jz) 

p_ = _ 2(m 2 + 7) 

2 >on P2on M{l + z) 

(i v. \— 2(/i 2 + 7 + 7t -2c ^ y^) , Ayl x 

(fc " ki) ™ = (A4) 
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where cosd = kj_ ■ kjj_, since we have chosen x = kjj_, 7 = k\ and 7; = k 2 ± . It should be 
pointed out that 1 > \z\ and 1 > \zi\ since we are considering only particles, i.e (p + /2)±k + > 
and (p + /2) ± kf > 0. 

In order to perform the analytic integration, one has to consider the following three poles 

fc __2(m 2 + 7) M . 2e 
L ~ M(l - z) ~Y~ 1 M(l-z) 
2(m 2 + 7 ) M . 2e 



ku M(l + z) + 2 + Z M(1 + z) 

2(7 + 7i + V 2 ~ 2cosd \/tTi) ■ 2e 



Ku ~ K M{z- Zi ) ^ u M{z- Zi )- 

M 2(7 + 7t + /i 2 - IcosQ v ^t") 2e 
~ 2 ^ M(z-Zi) + *M(z-^) { } 

First of all, let us mention that the cases \z\ = 1 and \zi\ = 1 lead to a vanishing 
contributions, as one can check by properly recombining the factors 1 / (1— z 2 ) and 1 / (z—Zi) in 
the integral ( 1A3I) and the previous poles. When 1 > z > Z{ > — 1 one can close the integration 
contour in the lower plane, taking the residue at k~ — k L - While for 1 > Zi > z > — 1 one 
can choose k~ = ky. Then, one has 

T {L) {l,z-^ h z u cos6) = ig 2 [ ^— — — 1— , x 



2tt M 2 (1-z 2 ) (z-Zi) 



1 1 



(Ar - k~ L ) (k~ - ky) {k~ - k 



LU) 



. 2 2 8 J 6(z- z, 

» 9 



M 2 (l-z 2 ) (z-Zi) 2 (m.2+ 7 ) — — 4- 2{m2+ ^ - M - j e 

V ' K % > V M(l-z) 2 ^ M(l+a) 2 ifc 
1 

2(m 2 + 7 ) _ M/i _i „ \ i_ 2(7+7j+M 2 ^2cos6> y/ry-) _ . 
M{\-z) 2^ + Zi ' + M(z-2i) ifc 

0(^-2) 



2(m 2 + 7 ) , M _ 2(m 2 +7) , M, • 
A/(l+z) 2 Af(l-z) 2 Lt 



1 



2(m 2 + 7 ) 1 M/i _ \ , 2(7+7z+At 2 -2cose y / 77l) . • 
' M(l+z) 2 ^ Z V M(z-Zi) Lt 



(A6) 



12 



Finally, one can rewrite 

X (L) (7, z\ 7*, z h cosO) = -g 2 



8 



1 



M 2 (l - z 2 ) (z - Zi ) M- 2(m2+7) - 2(m2+7) + ie 

K ' K ' m M(l+z) M{l-z) ^ lt 



X 



9{z - Zi ) 



M(l-z) 2 



f + * 



2(7+7 t +M 2 -2cos6> 
M(z-z;) 



2(m 2 +7) _|_ M-(\ — z -) + 2 (7+7i+<" 2 -2cos6> VT^) 



Af(l+z) 2 



M{z~Zi) 
1 



7 + (1 — z 2 )k 2 + z 2 m 2 — ie (z — Zj 
0(z - Zi) 



Q 2 +7) _ M£ 
(1-z) 4 



(z-Zi) 



ie 



.Zi - z 



(to 2 +7) _|_ M % (-\ _ „ \ 1 (7+7i+M 2 -2cos6» y/TT?) 
(1+z) + 4 



1 



+ ie 



7 + (1 — z 2 )k 2 + z 2 m 2 — ie 
where £^(7, 2; 7«, cosO) is given by 

<7 (L) (7, z\ 7i,Zi,cos9) 



(z-Zi) 

,z;~fi,Zi,cosV) 



(z - Zi) I ( m2 ±l) _ Ml(\ 4. r.) + (7+7»+At 2 -2cos6> y/^yT) _ 
I (1— 2) 4 \ */ (z— z.1 



(*—*») 



ie 



[Zi— z 



(to 2 +7) . M 2 ft _ \ , (7+7i+M 2 -2cos6» v^pyl) 
(1+z) 4 V 1 Z V 



0(z -*)(!-*) 



( z _ z .) [ m 2 + 7 _ + z .)(l _ z )] + (1 _ z )( 7 + 7 . + ^2 _ 2co ^ _ le 

6( Zi -z)(l + z) 



(z - Zi) [M-(l - + z) - (m 2 + 7)] + (1 + z)( 7 + % + /i 2 - 2cos9 ^tu) 



ie 



{z- Zi )(l-z) 



+ 



/3(z, z») + 7(1 - - 2(1 - z)cos0 v /tt~ - ie 

6(zj-z)(l + z) 

(3{-z, -Zi) + 7(1 + Zi) - 2(1 + z)cos8 y/yfi - ie 

where (3(z, Zj)is given by: 



(3(z,Zi) = (l-z) 



M 2 



ji + jj 2 - (z - Zi)(l + z^— m 2 



+ (1 - Zl )m 2 



(1-z) 



M 2 

/i 2 + — (l-z)(l + z l )-2m 2 



+ (1 - Zi)m 2 
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where M 2 /4 = (m 2 + 7$) /(l — z 2 ) has been used. For further purposes, the denominator can 
written as 



P(z,Zi) +7(1 - Zi) 



(1 - z l )[ 1 + (1 - ^)/€ 2 + ^m 2 ] + p(z, z^ - (1 - 2i )[(l - z*)k 2 + z 2 m 2 ] 



(1 _ Z .) A + (1 - z) 
{I- Zi )A + {l- z)B 



M 2 

V 2 + — (1 ~ z)(l + Zi ) - 2m 2 



+ {l-Zi){l-z 2 ) 



M 2 



(A10) 



where A = 7 + (1 - z 2 )k 2 + z 2 m 2 and B = /j 2 - 2k 2 - zziM 2 /2 

By using the Feynman trick to reduce X^ L ^(7, z; 7$, z*, cos#) to a single denominator, one 
has 



J (L) (7,2; 7^,21;, cos0) = g 2 da x 

■/ 

^-^(l-z) 



[A* + ((1 - Zi)A + (1 - z)B - 2(1 - z)cos6 ^77") (1 — a) — ie\ ' 

6{zi-z){l + z) 



+ 



[Aa + ((1 + Zi)A + (1 + z)B - 2(1 + z)cos# ^77") (1 
da 9(z — Zi)(l — z) 



a — le 



+9' 



o [1-^(1 -a)] 2 
da 



(l-z)(l-a) 



0(^-z)(l + z) 



+ 



o [1 + ^(1 -a)] 1 



2 ■ 



(All) 



where [1 ± Zj(l — a)] > since 1 > |zj| (|zj| = 1 can be excluded for physical motivation, 
since this amounts to an infinite energy for one of the incoming particle). The change of 
variable in the first integral is 



y 



and in the second one is 



(i-z)(i-q) 

[1 - z { (l - a)] 

(l + z)(l-q) 
[1 + Zi(l -a)} 
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Then, one has 



[ z - z .) e[l-z-y{l-Zi)] 



9" I dy 

>o \A + yB — 2ycos9 ^/yyl — ie] ' 



+ 



+ 9 2 f dy d ( z *- z ) d ^ + z-y^ + ^)] 

Jo 



[A + yB — 2ycos6 ^J 7 ^ — ie] 

1 [ [A + yB — 2ycos9 ^/tt7 — ie] 

| 6(z l -z)6[l + z-y(l + z l )} 
[A + yB — 2ycos8 v /tt7 — ie] 

Let us change y — » — y, then one gets 

7-i [ [A — + 2ycos6 y/yfi — ie] 

6{ Zl -z) 9[l + z + y(l + Zi)\ ' 



[A — yB + 2ycos6 y/yyl ~ & e ] ' 



dy / drf 



o(-y) Ki-iM) 



(A12) 



x 



-i J-oo [7' + 7 + (1 - z 2 )k 2 + z 2 m 2 + y^-zzi + 2ycos9 y/jjl - ie] ' 
{9{z - Zi )9[l-z + y(l - Zi)} + 6{Zi - z) 9 [1 + z + y(l + *)]} , (A13) 

where the term y(M 2 /2)zzi has been inserted for obtaining the same denominator in Eq. 
(159]) and 



M 2 

la{y) = ~yB - y—zzi = y(2n 2 - /i 2 ) 



(A14) 



Appendix B: The kernel V. in ladder approximation 



In the ladder approximation of the 3D LF equation for bound states (cf Eq. (150]) ) it is 



present the following kernel V t 

V b iL \ 7 ,z- n ',z')= i 2 g 2 p 



(L) 



dk~ 



d A k" 



1 



2vr [(E + k) 2 - m 2 + ie] [(§ - k) 2 - m 2 + ie] 
1 



x 



(2vr) 4 ^" 2 +p-k"z'-j'-K 



g 2 p' ] 



d A k" 



2 + le f (k-k») 2 -^ 2 + 
1 



ie 



(2VT) 4 [jfe/12 + p . _ y _ K 2 + ie ] 3 

1 1 



dk' 



27T [(§ + *) 



2 _ m 2 



ie] [(I -£;) 2 -m 2 + ie] {k - k") 2 - fi 2 + 



(Bl) 
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An explicit expression can be found in Ref. Another explicit expression can 

be obtained from as given in subsect. IV C\ and recalling that V&(7, z\ 7', CO = 

lim c ^o K(7, z] li, Zi, 7', C, C, cos9). 

In this Appendix, it will be illustrated in detail the kernel for the scattering-state equation 
(cf Eq. (I69p ). In this case, the kernel acquaints new dependencies: i) upon another compact 
variable z" and ii) the cosO = • kjx- Moreover, there is an understood dependence upon 
7i, Zi and p. In ladder approximation, one has the following expression for V S {L) 

"°° dk~ 



-p' 



2tt 



Gf\k,p) x 



d A k" 



& L \k,k",p) 



(2k) 



k" 2 + ^ -m 2 + p- k" C + 2k" -kiC-l' + ie 



d 4 k" 



dk~ 1 1 

~2t7 [(2 + k) 2 -m 2 + ie] [(§ - k) 2 - m 2 + ie] 
1 1 



(2vr) 4 [A;" 2 + p ■ fc"C' + 2A;" • fcjC - 7' - 
dfc" 1 



0V 



2 + ze ] 3 (* - A;") 2 - /i 2 + 
1 



2n [(§ + k) 2 - m 2 + ie] [(§ - fc) 



2 _ m 2 



(B2) 



where 



p(fc, 7 / ,c,c / ,co^)= y 



(fc - A;") 2 - /i 2 + 



Adopting the formula (see also Ref. 



3v 2 dv 



ab 3 J [a(l - u) + 6u] ■ 



where 



■ l3 



(2tt) 4 [p 2 + p . k"(' + 2k" ■ hC - i - k 2 + ie] 



x 



(B3) 



(B4) 



a — (k — k"y - p 2 + ie 



b = k" 2 + p ■ k"(' + 2k" ■ ki( - 7 - k 2 + ie 



(1 - v)a + vb={l- v)[{k - k") 2 - fi 2 ] + v[k" 2 + p ■ k"C + 2k" ■ k { ( - 7' - k 2 } 



k" 2 - 2k" ■ [(1 - v)k - <| - vkiC] + (1 - v){k 2 - p 2 ) - v(i + k 2 ) + ze 
q 2 - [(1 - - <| - ^C] 2 + (1 - v)(k 2 - p 2 ) - v{i + k 2 ) + 



(B5) 
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with 



Then one has 



q — k" — (1 — v)k + v('- + vkiC 



V(k,j',(,(',cost 
-i 



(2vr) 4 J 
v 2 dv 



3v 2 dv 



d 4 q 



[q> + A(kXA',cose)+ie} 



2(4vr) 2 y [A(k,(,(',cos9) + ie] 



(B6) 



where 



A(k,CX',cosi 



fc 2 v(l - v) + 2(1 -v)v k ■ [C'| + fcCl - ^[C'f + **C] 2 - (1 - ^)/i 2 - «(Y + «: 2 ) 



f (1 — v) 



k~ [k + + C'— - ^C^r 



M. 



z( -zC < —k i -i-2(,cos0y/Tfi 



-v 2 [('l + k l (f-(l-v)^-v( 1 ' + K 2 ) = 
= v(l - v)k~ (k+ + C'y - ZiC^pj - £ D (v, z, T ', C, (', cost 
with (dependence upon 7^ z t and are understood) 



£ D (v, z, 7', C, C', cos6) = v(l - v) 



M 2 M 

—z( +zC,—k i + 7 + 2Ccos^v^77 



(B7) 



+ ^ 2 [C'f + ^C] 2 + 



+(1 - v)fi 2 + u(V + K 2 ) 
"M 2 



f (1 — 



X' + 7 



+ ^ 2 (C ,2 ^ + ^C 2 ) + 



+ (1 - w)/i 2 + v(i + K 2 ) + V{\ - v)( 
' M 2 



' M 2 „ 
■^i - ^ — h 2cos6 , A /77i 



u(l — u) 



-zC' + 7 + « 



C' 2 ^ + « 2 (i + C 2 ) 



+ 



+ (l-w)/i 2 + ^7' + w(l-i;)C 



M 2 „ , 

zz %—^~ + 2cos6'- v /777 



(B8) 



Let us recall that p ■ ki = leads to /c/" + fc^ =0 and 7; = M 2 (l — -2 2 )/4 — m 2 as mentioned 
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in the previous Appendix. Then the kernel for obtaining the scattering state becomes 



vj L \ 1 ,z ]li ,z i ^',c,C',cose) = 

2 + % r°° dk~ 1 1 

~ 9 P 2(4tt) 2 J /_ 00 ~27 [(| + A;) 2 -m 2 + ze] [(| - k) 2 - m 2 + ie] * 

w 2 rfw 2 . i 1 r 1 2 , f° dfc" 

5 P nTTjo , a / . , : x , m — r~r v dv — x 



[A(p, fc, C, C, cos0) + ie} 2 y ' 2(4tt) 2 (f + k+) (f - fc+) J 2vr 

1 1 

x 



[(I + k)~ - (| + k)~ n + ie/(M/2 + fc+)] [(I - fc)~ - (| - fc)- + ie/(M/2 - fc+)] 

1 

[«(1 - v)k- (k+ + C'f - *Ct) - 7, 7', C, C co^) + ie] 2 ~ 

2 % 1 f 1 2 r f°° dk ~ 

= - g —r —r- r— / v dv — — x 



2tt 2 M 3 (l - 2 2 ) i ./_«, 2tt 

1 



[(I + k)~ - (§ + fc)~ + <2e/[M(l - *)] [(1 - - (1 - + z2e/[M(l + *)])] 

1 

{v(l -v)((>-z- ZiO [k- -k D + i2e/[Mv(l - v)((> - z - z^)}] } 2 ~ 
where 



(B9) 



(P + b)~ = ~ 2 ( m2 + 7) 
l 2 >on Plon M{l-z) 

(t-lc\- = - 2 ( m2 + l) 

2 )on P2on M{l + z) 
2 £ D (v,7,z,7',(,C,cos9) 

k D = ~TTT\ \T7< TV B1 ° 

Mv(l - v)(C - z - ZiQ 



One has the following three poles (1 > z > — 1) 



fa = t °- V-,'-^) (B11) 
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If (' > z + Zi( one can close in the upper plane taking the residue at k p , i.e. 



1 



1 



dk 

~^ [(§ + W - (I + k)on + ie] [(I - k)~ - (| - k)~ n + iej] 
1 



x 



[v(l-v)((> -z-z^ik- -k D ) + ie\ 
1 



—i 



X 



1 



{«(1 - w)(C - * - *C)[f - (f - *ta - k D ] + ze} 2 
M 1 



v(l - v)((> - z - Zi ()(f -±a*g.-kB) + 



\ 2 m 1+2 

If z + Zi( > £' one can close in the lower plane taking the residue at k n , i.e. 
dk~ 1 1 



(B12) 



2tt [(§ + *)" -(§ + *)" +ie] [(§_*)" 
1 



1 



X 



[ v (l -A£-ie)]' 

1 

Hi - tO(c - i - icFf + (I + *0™ - *b - ^]} 2 

. M 1 



[M-(f + -(§-*)- + fe] 



x 



(1-2 2 ) 



av + ie 



vil-vXC-z-ztO ( 



M 2 m 2 +7 

2 A/ (1-2) 



+ A; D + ie) 



(B13) 



Then, by exploiting Eq. (IBlOp one has 

Vj L \ 1 ,z- li ,z l , 1 ',C,C , ,cos9) = 

/ dv v 2 F(v, 7, z; 7', C, C, cos6» 
Jo 



2(2vr) 2 (1 - 2 2 



x 



2(4vr) 2 [7 + (1 - z 2 )k 2 + z 2 m 2 - ie] J 



dvv 2 F{v, 1 ,z ] i,C,C\cos6) (B14) 
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where 



F(v, 7, z; Y , C, C', cost 



9((' - z - ztC) 



v(l - v)(C - z - ZiO - - £ D (v, 7, z, 7', C, C cosO) + *e 

0(z + ZiC - C) 



4 



cosy + 



+ 



;(1 -v)(z+ Zi ( - CO (¥ - - 7, *, y, c, e, 

(1 + z) 2 fl(C - Z - ZjQ 

[D(v, z, z u C, CO + v(l -v)(l + z)( (^zzi + 2cos6^rn) - ie] ' 

(1 - zf 6(z + ~ CO 

[V{v, -z, - Zi , C, -CO + v(l -v){l- z)C (^zzi + 2cos9^al) - ie] ' 



with 



V(v, z, z h C, CO = v(l - v)(C - z - ZiC) [m 2 + 7 - (1 + z) 



M 7 



+ 



f M 2 

+ (1 + z)£ D {y^ } z^\CX\cos6) -v(l - v)(l + z)( \— zz i + 2cos6' v /t7 



= v(l — v) 
+ (1 + z) ji; 2 
= v(l — v) 



,M 



(C'-z- Zi C)(k 2 + 7) + z 2 — (1 + C' - CZi) + (l + z)(j + K ) 



C' 2 ^ + - 2 (i + C 2 ) 



4 

+ (1 - v)[i 2 + U7 ; 
M 2 



(1 + C' - z t C)(n 2 + 7) + — z 2 (l +C- Czi) 



+ 



+ (1 + z) <N 



C' 2 ^ + * 2 (i + C 2 ) 



+ (i - v)^ 2 + «y 



X(v, Zi , C, CO [7 + z 2 m 2 + k 2 (1- z 2 ) + Y{v, z, Zi , C, C', 70] 



In Eq. ( 1B16j) . the following notation has been used 



Notice that 



x(v, = v(i-v)(i + {'-*{) 

T(v,z,z t ,C,C',l , )= .. (1 + Z) x 



(1-v) 



(1 + C - ZiC) 

,,2M 2 2 2 , 

C ^ + ^(l + C) + 7 



H h 7 



->x(^,c,C',)>o 



->x(u,-^c,-C',)>o 



where the positivity follows from the presence of the theta functions in Eq. (IB15j) 
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Then one can rewrite 



F(v, 7, z; i , C, (', cost 



(1 + zf 



x 2 (v,z h C(>) 

9(C - z - z t () 



■ l2 



[ 7 + Z 2 m 2 + K 2 (1 _ Z 2) + r(Uj Zj ^ C C / ; y) + Z ^ C> C /. ^ _ i£ ] 
^X 2 (v,-z h (,-(>) 

e{z + ZiC-C) 



where 



and 



[ 7 + Z 2 m 2 + K 2 (1 _ 2 2) + r(Uj _ Zj _^ C) _ C / 5 y) + Cj _ (f . _ Z . )C _ le] 



M 2 

C = — z Zi + 2cos6y/yy~i 



(B19) 



z(z,(,C;z l ) = ( 



(l + z) 



(B20) 



(l + C'-^C) 

Once more, the presence of the theta functions in Eq. (IB15I) helps to find the following 
constrain 



|z(*,C,C',*)l<ICI 



\z(-z,c,-C',-Zi)\ < ICI 



(B21) 



Notice that one can always put e — > X(v, Zi, (, (') e, given the limits on X(v, Zi, (, (')■ 

The kernel vj L ^ can be evaluated by slightly elaborating the standard Feynman trick. 
Indeed one can use 



1 



lim — 



B A 2 a^o A [BA B(A + X)_ 



lim — 

A^O A 



1 



1 



[B - i{B - A)} 2 Jo ^ [B - i{B - A) + £A] 2 



(B22) 



By applying the above relation with the following identifications 

At = 7 + z 2 m 2 + k 2 (1 - z 2 ) + T(v, z, z h C, (', 7) + Z(z, C, C'l Zi)C - te 

B± = 7 + z 2 m 2 + k 2 (1 - z 2 ) - te (B23) 



one gets 



(l + z) 



x 



l(l + (>- Zi () 
(1-*) 

+ (i-C + ^C) 



2(4tt) 2 

9(C - z - ZiC) H'(j, z, zf, y, C, (', cos6, /i 2 )+ 



0^ - c' + ZiQ -z, - Zi ] y , c, -c', ^ 



(B24) 
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where 



H'(j, z, Zi\ 7', C, (', cos9, /i 2 ) = 

= lim i [ ft( 7 , z, zf, i , C, C, co^, /i 2 , A) - ft( 7 , 2:, 7', C, C, cos0, /x 2 , 0)] (B25) 



with 



tt( 7) z, Zi-, V, C, C, coag, /i 2 , A) = ^ (1 + ^ ^ I TT^To f\ 



(l + S'-ZiQJo (l-v) 2 Jo 

1 



{7 + 2 2 m 2 + K 2 (1 _ ^ + £ [ r (^ Zj ^ Cj C) y) + ^ C) C . Z . )C + A] - se}' 

d£ / di' x 



(l + C-^Oio (l-^) 5 JO 
[7 + z 2 m 2 + k 2 (1 - z 2 ) + 7" + f C, C; - ie] 



(B26) 



It should be pointed out that, in order to obtain Eq. (1B24[) . the limit on A has been 
exchanged with the integral over v. 

By inserting the following change of variable 



z' = zz(z,(,C;zi 



(B27) 



such that \z'\ < |C| < 1 from Eq. (1B2 1|) . one gets 
H(y, z, Zi] 7', C, C', cos8, /1 2 , A) - 



dv 



{1 + C-ZiQJo {1-vY 
e(z') 6(z(z, C, C; *) - z') - 6{-z') d{z> - z{z, C, C; *)) 



d 7 " 



X 



[7 + z 2 m 2 + k 2 (1 - z 2 ) + 7" + z'C - ie 
(1 + z) f l dv 



5 



' „ ,r(v,z,z i ,t,t'>'f) + A 
7 — z - 



z(z,c,C';^ 



Z 1 jz(*,c,c';*)l x 
rf7 Z./ 2 z(*,c,C';*) 



(l + C'-^Oio (i-^) 2 , 
5 [ 7 "z(z, C, C; *) - z'T(v, z , z u c, C, i) - z'\] x 
0(z') e(z{z, c, C; *) - z') - e(-z') e(z> - z(z, c, C; *)) 



[7 + z 2 m 2 + k 2 (1 - z 2 ) + 7" + z'C - ief 



(B28) 
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Finally, let us perform the change v — > y/1 + y, then one has 

(1 + *) 



H(-f,z,Zi;j ,C,t,cos6,fi , A) 
•i 



oo roo 



dz' 5 



(i + C'-*C) 
V'^fo C, C'; *») - z'f( y , z, c, C', V) - z'x 



dy I d^" x 



x 



e(z') 6(z(z, c, C; - + 0(-zO ^ - z(z, c, C; *)) 



• l2 



[7 + z 2 m 2 + k 2 (1 - z 2 ) + 7" + z'C - ze] 

d // Z" 1 ^ g(fO C, C ; ; ti) - *') + g(-*0 0(z' - c, C ; ; *)) 

[7 + z 2 m 2 + k 2 (1 - z 2 ) + 7" + z'C - ief 



00 j — 1 



dy 5 
d 7 " 



7 C-^r(?/,c,C,7) — - — ; — a 



(1+*) 

1 d^ 0^) 0(z(z, c, C; *) - «') - 0(-*O 0(z' - z(z, c, C; *)) 
1 



• t2 



[7 + z 2 m 2 + k 2 (1 - z 2 ) + 7" + z'C - ze] 



-(y*A + yB + ^) 

y 



where 



r(y,(,C,i)= y 



C' 2 ^ + - 2 (i + C 2 ) 



+ ^^ + (l + y)V 



and the argument of the delta has been recast in the following form 

r(y,C,C,7) + — 77—^ — A -7 - = 
(1 + z) z' 

= - [y 2 A{C C, 7', + 2/ B(z, z', z h C, C, V, 7", /A A) + /i 2 ] 



with 



^(C,C / ,7 , ^ 2 ) = C ,2 ^ + ^ 2 (l + C 2 ) + 7 , 

z', zt, C, C, 7', 7", /A A) = /i 2 + i - + ^ 



M 2 



(1 + *) 



The integral over y is readily done, obtaining 
"H(7, z, 2^; 7', C, cos6, /i 2 , A) 

1 + 2 



d 7 " 



1 dz' 



X 







1 z ' 

0(z') (9(C) - 9(-z') 6{-Q 



1 + C-ZiC C y [7 + z 2 m 2 + k 2 (1 - z 2 ) + 7" + z'C - ze] 2 



ft 2 ) — /l 2 j 
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where we have used Eq. (IB20j) (exploiting the theta functions for determining the signs of 
C) an d 

with /(^j) = and the derivative evaluated at y = y^. Notice that 6(yi(X)) comes 

from the lower extrema in the integral on y. Moreover, for simplifying the notation, we have 
dropped the dependence upon some variables in y±, that represent the two solutions of the 
second order equation 

y 2 A(C, C, 7', k 2 ) + y B(z, z', z h (, (', 7', 7", /A X) + fj 2 = (B35) 

Once "H(7, z, z\\ 7', C, C'> cos6*, fi 2 , X) is formally obtained, its derivative with respect to A, at 
A = 0, can be evaluated as follows 

™oo /•! AJ 

X 



/ / / f°° f dz' 

'H'(l,z,z i ;-f',C,C',cos9,fi 2 ) = dr/'j — 

J-oc J-l Z 

1 + z z'\ e(z') 9(C) ~0{-z') 6(-() 



1 + C - ZiC C / [7 + z 2 m 2 + k 2 {1 - z 2 ) + 7" + z'C - ie] 



2 



k 2 ) — ^ 2 | ^ 



-%(0)) 2/i 2 yi(0) 



(2/ 2 (o)^(c, c, y, * 2 ) - ^WMM, e, y, « 2 ) - » 2 \ 

00 .7" /" dz> ^^W', Z 'XX,cosQ^) 



'DC: 



1 [7 + z 2 m 2 + K 2 (l-z 2 )+7" + ^C-ze] 



where 



with 



Q > ( z, z, y , 7", y, c, c, , 2 ) - * (t^^ - f) ^ m ^ 9{ -° x 

a(z,^,C,CW,7V*,/) (B37) 



ym 



~ %l(0)) (y 2 (oM(c,c',y,^ 2 ) ^MSCcv?^ -v 2 \] (B38) 
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and 



y±(o) 



1 

2^l(C,C',7 / ,^ 



x 



-B{z, z\ C, C, Y , 7", /A 0) ± V, V h C, C, i , 7", j* 2 , o) - V -4(C, C, Y , 

£*<0) - T (1 + C/ -f° ^ (0) (B39) 

^ (1 + *) v/i3 2 (z, z', ^ C, C, Y, 7", ^ 2 , 0) - WW, C, Y, 

Let us recall that Ui(0) still depends upon z,Zi,( and and therefore the proper changes, 
dictated by the arguments of Q', must be performed in order to get the final expression for 
Vs . The necessity of real solutions for getting non vanishing Nakanishi amplitude imposes 
a constraint on the discriminant of Eq. (1B35[) , viz 



b\z, z', Zi , c, C, Y , 7", /A o) - V AC, C, Y , 



H + 7 - 7 — j -4/i 



C' 2 ^ + - 2 (i + C 2 ) + Y 



> 



(B40) 



More constraints can be obtained from the positivity of the solutions Vi(0). In particular 
one has 

• if -4(C, C') 7'; > 0, then one has to consider B(z, z', Zi, C, C'; 7'; 7") ^ 2 ; 0) < 

7', k 2 ). In this case, two solution for ^(0) are allowed; 

• if A((, C, 7', k 2 ) < 0, then B(z, z\ z i} (, 7', 7", /x 2 , 0) > 0. Only one solution is 
allowed. 



In conclusion, from Eq. ( 1B24|) and using Eq. ( 1B36[) . one has 



K (L) (7,^7 l ,^,7",C,C', 



cost 



2(4?r)^ 



drf j dz' X 



1 



[7 + z 2 m? + k 2 (1 - z 2 ) +<y' + z' {^zzi + 2cos#^77~) - ie] 
(1 + z) 



L(i + c-ZiC) 

(l-z) 

+ 



0(C - z - *C) 7", Y, Y, c, C, /i 2 )+ 
z - C' + ZiC) Q'(-z, 7", Y, z'C, -C, ^) 



(B41) 



;i-c+*o 

where for the sake of notation adopted in Sect. HVI the change 7" — > 7' has been performed. 



[1] E.E. Salpeter, H.A. Bethe, Phys. Rev. 84, 1232 (1951). 



55 



[2] N. Nakanishi, Prog. Theor. Phys. Suppl. 43, 1 (1969); 95, 1 (1988). 
[3] G.C. Wick, Phys. Rev. 96, 1124 (1954). 

[4] S.M. Dorkin, L.P. Kaptari, C. Ciofi degli Atti, B. Kampfer, Few Body Syst. 49, 233 (2011). 

[5] F. Gross and A. Stadler, Few Body Syst. 49, 91 (2011). 

[6] C. D. Roberts, Prog. Part. Nucl. Phys. 61, 50 (2008). 

[7] K Nishijima, Prog. Theor. Phys. 14, 203 (1955). 

[8] C. Schwartz and C. Zemach, Phys. Rev. 141, 1454 (1966). 

[9] S. G. Bondarenko, V. V. Burov, M. Beyer and S. M. Dorkin, Few-Body Syst. 26, 185 (1999); 

S. M. Dorkin, M. Beyer, S. S. Semikh and L. P. Kaptari, Few-Body Syst. 42, 1 (2008). 
[10] N. Nakanishi, Phys. Rev. 127, 1380 (1962); Phys. Rev. 130, 1230 (1963) Graph Theory and 

Feynman Integrals (Gordon and Breach, New York, 1971). 
[11] V.A. Karmanov and J. Carbonell, Eur. Phys. J. A27, 1 (2006). 
[12] J. Carbonell , V.A. Karmanov, Eur. Phys. J. A27, 11 (2006). 

[13] K. Kusaka, A.G. Williams, Phys. Rev. D 51, 7026 (1995); K. Kusaka, K. Simpson, 

A.G. Williams, Phys. Rev. D 56, 5071 (1997). 
[14] N. Nakanishi,Phys. Rev. 133, B214 (1964); Phys. Rev. 135, B1224 (1964). 
[15] S. J. Brodsky, C.-R. Ji and M. Sawicki, Phys. Rev. D 32, 1530 (1985). 
[16] Dae Sung Hwang and V. A. Karmanov, Nucl. Phys. 696, 413 (2004). 

[17] V. Sauli, J. Adam , Phys. Rev.D 67, 085007 (2003); V. Sauli, J. Phys. G 35, 035005 (2008) 

[18] J. Carbonell, V.A. Karmanov, M. Mangin-Brinet, Eur. Phys. J. A39, 53 (2009). 

[19] J. Carbonell, V.A. Karmanov Eur. Phys. J. A46, (2010) 387. 

[20] J. Carbonell and V. A. Karmanov, Few-body Syst. 49, 205 (2011). 

[21] P. C. Magalhaes, M. R. Robilotta, K. S. F. F. Guimaraes, T. Frederico, W. de Paula, I. Bedi- 
aga, A. C. dos Reis, "Towards three-body unitarity in D + — > K~tt + tt + '\ arXiv:1105.5120v2 
[hep-ph], in press Phys. Rev. D. 

[22] G. F. Giudice, R. Rattazzi, J. D. Wells, Nuc. Phys. B 630, 293 (2002). 

[23] J. Carbonell, B. Desplanques, V.A. Karmanov and J.F. Mathiot, Phys. Reports, 300, 215 
(1998). 

[24] J. H. O. Sales, T. Frederico, B. V. Carlson, P. U. Sauer, Phys. Rev. C 61, 044003 (2000). 
[25] J. H. O. Sales, T. Frederico, B. V. Carlson, P. U. Sauer, Phys. Rev. C 63, 064003 (2001). 
[26] T. Frederico, J. H. O. Sales, B.V. Carlson and P.U. Sauer Nucl. Phys. A 737, 260c (2004). 

56 



[27] J. A. O. Marinho, T. Frederico, P. U. Sauer, Phys. Rev. D 76, 096001 (2007). 
[28] J.A.O. Marinho, T. Frederico, E. Pace, G. Salme, P. U. Sauer, Phys. Rev. D 77, 116010 
(2008). 

[29] T. Frederico and G. Salme, Few-body Syst. 49, 163 (2011). 
[30] R.E. Cutkosky, Phys. Rev. 96, 1135 (1954). 

[31] S. J. Brodsky, H. C. Pauli and S. S. Pinsky, Phys. Rep. 301, 299 (1998) 
|arXiv:hep-ph/9705477] . 

[32] S. J. Brodsky, C. D. Roberts, R. Shrock and P. C. Tandy, Phys. Rev. C 82, 022201 (2010). 

[33] P. A. M. Dirac, Rev. Mod. Phys. 21, 392 (1949). 

[34] R. M. Woloshyn and A. D. Jackson, Nucl. Phys. B 64, 269 (1973). 

[35] C. Itzykson, J.B. Zuber "Quantum Field Theory", Dover Publications (2006). 

[36] B.L.G. Bakker, J.K. Boomsma and C.-R. Ji, Phys. Rev. D 75, 065010 (2007). 

[37] J.R. Cooke, G.A. Miller and D.R. Phillips, Phys. Rev. C 61, 064005 (2000). 

[38] M. Mangin-Brinet and J. Carbonell, Phys. Lett. B 474, 237 (2000). 

[39] F. Gross " Relativistic Quantum Mechanics and Field Theory", John Wiley and Sons, New- 
York (1993). 

[40] N.E. Ligterink and B.L.G. Bakker, Phys. Rev. D 52, 5954 (1995). 
[41] M. Sawicki, Phys. Rev. D 32, 2666 (1985); D 33, 1103 (1986). 

[42] A. Castro Neto, F. Guinea, N. M. R.Peres, K. Novoselov, and A. Geim, Rev. Mod. Phys. 81, 
109 (2009). 

[43] N. M. R. Peres, Rev. Mod. Phys. 82, 2673 (2010). 



57 



